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PREFACE 


I T is now over thirty years since Heaviside’s operational methods of 
solving the differential equations of physics were first published, but 
hitherto they have received very little attention from mathematical 
physicists in general. The chief reason for this lies, I think, in the lack 
of a connected account of the methods. Heaviside’s own work is not 
systematically arranged, and in places its meaning is not very clear. 
Bromwich’s discussion of his method by means of the theory of functions 
of a complex variable established its validity; and as a matter of practical 
convenience there can be little doubt that the operational method is far 
the best for dealing with the class of problems concerned. It is often said 
that it will solve no problem that cannot be solved otherwise. Whether 
this is true would be difficult to say ; but it is certain that in a very 
large class of cases the operational method will give the answer in a 
page when ordinary methods take five pages, and also that it gives the 
correct answer when the ordinary methods, through human fallibility, 
are liable to give a wrong one. In particular, when we discuss the small 
oscillations of a dynamical system with n degrees of freedom by the 
method of normal coordinates, we obtain a determinantal equation of 
the wth degree to give the speeds of the normal modes. To find the ratios 
of the amplitudes we must then complete the solution for each mode. If 
we want the actual motion due to a given initial disturbance we must 
solve a further family of 2n simultaneous equations, unless special 
simplifying circumstances are present. In the operational method a 
formal operational solution is obtained with the same amount of trouble 
as is needed to give the period equation in the ordinary method, and 
from this the complete solution is obtainable at once by a general rule 
of interpretation. For continuous systems the advantage of the opera- 
tional method is even greater, for it gives both periods and amplitudes 
easily in problems where the amplitudes cannot be found by the ordinary 
method without a knowledge of some theorem of expansion in normal 
functions analogous to Fourier’s theorem. Heat conduction is also 
especially conveniently treated by operational methods. 

Since Bromwich’s discussion it has often been said that the operational 
method is only a shorthand way of writing contour integrals. It may be ; 
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but at least one may reply that a shorthand that avoids the necessity 

of writing — J dK in every line of the work is worth while. Connected 

with the saving of writing, and perhaps largely because of it, is the fact 
that the operational mode of attack seems much the more natural when 
one has any familiarity with it. After all, the use of contour integrals 
in this connexion was introduced by Bromwich, who has repeatedly de- 
clared that the direct operational method of solution is the better of 
the two. 

My own reason for writing the present work is mainly that I have 
found Heaviside’s methods useful in papers already published, and shall 
probably do so again soon, and think that an accessible account of them 
may be equally useful to others. In one respect I must offer an apology 
to the reader. Heaviside developed his methods mostly in relation to the 
theory of electromagnetic waves. Having myself no qualification to write 
about electromagnetic waves 1 have refrained from doing so; but as the 
operators occurring in the theory of these waves are mostly of types 
treated here I think the loss will not be serious. It can in any case be 
remedied by reading Heaviside’s works or some of the papers in the list 
at the end of this tract. 

A chapter on dispersion has been included. The operational solution 
can be translated instantly into a complex integral adapted for evalua- 
tion by the method of steepest descents; a short account of the latter 
method has also been given, because it is not at present very accessible, 
and is often incorrectly believed to be more difficult than the method of 
stationary phase. Two cases where the Kelvin first approximation to the 
wave form breaks down are also discussed. 

My indebtedness to the writings of Dr Bromwich is evident from the 
references in the text. In addition, the problems of 4.4 and 4.5 are taken 
directly from his lecture notes, and several others are included largely 
as a result of conversations with him. 

My thanks are also due to the staff of the Cambridge University 
Press for their care and consideration during publication. 

HAROLD JEFFREYS 


St John’s College, 
Cambridge. 

1927 July 19. 
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CHAPTER I 

FUNDAMENTAL NOTIONS 


1 . 1 . Let us consider the linear differential equation of the first order 

« 

where R and S are known functions of bounded and integrable when 
O^x^a. Suppose further that ^ = when ^ = 0. Let Q denote the 
operation of integrating with regard to x from 0 to x, so that 

Qy^^^ydx (2) 

Perform the operation ^ on both sides of the equation (1). Then we 
find 

y-y,--Q{Ry^ S\ (3) 

and the right side vanishes with x. This can be rewritten in the forms 

{l-^qR)y^y,^QS, (4) 

QS+QRy (5) 

These are botli equivalent to the original differential equation (1) to- 
gether with the given terminal condition. When we write QRy we mean, 
of course, that R is to be multiplied into y and the product integrated 
with regard to x from 0 to x. But the whole expression for y may be 
substituted in the last term of (5), giving in succession 

y-^y.^QS^-QJt{y,^-QS+qRy) 

-yo^qS^qR (yo + qs) + qRqR ( 3/0 + qs-^ qRy) 

= (yo + qs) + qR ( 3/0 + qs) + qRqR {y, + qs) 

+ qRqRqR{yo + qs)^ ( 6 ) 

on repeating the substitution indefinitely. We have to show that the 
infinite series (6) converges and that it is the correct solution. In 
evaluating each term R is supposed to be multiplied into the whole 
expression after it, and q to operate on the whole expression after it. 
Suppose that within the range considered 

\R\^A, \yo+qS\^B, (7) 

where A and B are finite. Then the absolute value of the second term 
J I 
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is less than ABx^ that of the third than A^Bx^l2\, that of the fourth 
than A^Ba^lS ! , while the general term is less than A^Bx^jn ! . The series 
therefore converges at least as fast as the power series for exp Ax, It 
therefore represents a definite function, and on substituting it in (5) we 
see that the equation is satisfied for all values of x\ also the solution 
reduces to when x=^0^ as it should. Thus (6) is the correct solution. 

The solution can also be written 

y^il^QR^QRQR^QRQRQR^ + QS) (8) 

The operator between the first pair of brackets is the binomial expansion 
of (1 - QR)~^^ carried out as if QR was merely a number. Since y^ + QS 
is a determinate function, we can write the solution in the form 

y = r:r^(3'« + Q^. (9) 

provided that yo + QS is evaluated first and that the operator (1 - QR)"^ 
is expanded by the binomial theorem before interpretation. In fact (9) 
is merely a shorthand rule for writing (8). But on returning to (4) we 
see that (9) is also the solution of (4) carried out as if 1 ~ QR was a mere 
number. 

It is evident from (8) that the values of 8 for negative values of x do 
not affect the solution provided y^ is kept the same. Suppose then that 
8 was zero for all negative values of x, and that y was zero when x-0. 
The solution would be 

y^{i^qR^qRqR^,..)qs, (lo) 

and this solution would be unaltered if the lower limits of all the 
integrals were replaced by - co . But if now we add to ^ a constant y^j^ 
for all values of x between 0 and qS will be increased by yo for all 
values of x greater than and (10) will be converted into (8). If then ^ 
tends to zero, y^ remaining the same, y^ will tend to yo» and we recover 
the solution (8) with the original initial condition. The physical interest 
of this result is that it corresponds to our notions of causality. Suppose 
that the independent variable x is the time, and that y represents the 
departure of some variable from its equilibrium value ; then R represents 
a property of the system and 8 an external disturbing influence. If the 
system was originally in its equilibrium state, the form (10) exhibits 
the disturbance produced by the external influence after it enters. If it 
was undisturbed up to time zero, the part of (8) depending on 8 repre- 
sents the effect of the finite disturbances acting at subsequent times, 
while the part depending on y^ represents the effect of the impulsive 
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disturbance at time 0 required to change y suddenly from zero to If 
we like we can separate the solution into two parts 

{l-QBY^y, and (l-^QE^QS, 

and say that the first represents the effect of the initial conditions and 
the second that of the subsequent disturbances. 


1 . 11 . The method just given can be extended easily to cover many 
equations of the second and higher orders. Thus if our equation is 

(■) 

with and ^ = yi when a: = 0, we find by integration 


-yi = Q(Ry+S), 


y-y^-a;yi = QQ(Jty + S) (3) 

This leads to the solution in a series 

y=(i + Q‘Ji + +. . .) (yo + + W (4) 

where j^/(t)didt. (6) 


1.12. As a special example, consider the equation 

(■) 

with y = 1 when x = 0. Then carrying out the process of 1.1 (6) we get 

^ = (1 + aQ + ...) 1 (2) 


= l + cU?4--^ + ..., 


the ordinary expansion of exp <ur. Or take the equation 



with y=l and dyjdx = 0 when a? = 0. We can infer 






the ordinary expansion of Jo{x), 


1-2 
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1 . 2 . The foregoing method is due originally to J. Caqu6*; it is a 
valuable practical method of obtaining numerical solutions of linear 
differential equations. Its extension to equations of any order, or to 
families of simultaneous equations of the first order, is more difficult 
unless special simplifications enter. But if the equations have constant 
coefficients, which is an extremely common case in physical applications, 
a considerable development is possible. This arises from the fact that 
the operator Q obeys the fundamental laws of algebra. Thus if a is a 


constant and u and v known functions of 

Q{au)=^a{qu\ (1) 

Q{u^rv) = Qu+ Qv, (2) 

= Q^Q^u - ( 3 ) 

Consequently Q behaves in algebraic transformations just like a number. 
If for instance we have two operators/(Q) and ff(Q) both expressible 
as sums of integral powers of Q, thus ^ 

^(Q) = ao + aiQ + a.,Q^ + a.^(^-^ ' ( 4 ) 

9(Q)-^o + hQ’^ 1^20^ ( 5 ) 


where the a^s and b*s are constants, let us for a moment replace $ by a 
number z small enough to make the series converge absolutely. Form 


the product series 

f{z)g(z) = (ao + aiZ + a2Z^+...)(bo + biZ-¥b2Z^+...) 

^Co-^CiZ + C2Z^’^ (6) 

say. Consider first the case where the series are both polynomials. Then 
if S is an integrable function of x 

f{Q)g{Q)S-{c,^c,q-^c2q^-^-.,.)S, ( 7 ) 

For the left side means 

{(Iq + (351 q + 6^2 + •••) (b^s + bi qs + ^2 q^^ + . . .) 

=^ao{boS+b,qS+b 2 q'^S+,..) + a^{boqS + b,Q^S-^...) 
^a2(b,q^S (8) 

* LiouvilWs Journal, (2), 9 (1864), 185-222. Further developments are given by 


Fuchs, Ann, d. Matem., (2), 4 (1870), 36-49; Peano, Math. Ann., 32 (1888), 450-456; 
H. F. Baker, Proc. Lond. Math. Soc„ (1), 34 (1902), 347-360; (1), 35 (1902), 333- 
378; (2), 2 (1905), 293-296; Phil. Trans., A, 216 (1916), 129-186. Caqu^ considers 
only a single differential equation, but notices that the operator in the result is 
in the form of a binomial expansion. For other operational methods based on these 
principles, but applicable to equations of order higher than the first, or to families 
of equations of the first order, the above papers of Prof. Baker may be consulted. 
Physical applications are given by W. L. Cowley and H. Levy, Phil. Mag., 41 
(1921), 584-607; Jeffreys, Proc. Lond. Math. Soc., (2), 23 (1924), 454 and 465;, 
M.N.R.A.S., Geoph. Suppl. 1 (1926), 380-383. 
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using (1) and (2); and then using (3) we can collect the terms involving 


the same power of Q and obtain 

aoboS+ (aobi + ai h) QS 4- + «i bi + bo) Q^S-h ... (9) 

which is by definition the same as 

(co + CiQ + C2Q"+...)^. (10) 


When the series are infinite, their convergence may be established 
easily. If the series for / (z) converges absolutely when \z\^r^ then a 
number M must exist such that anr''^ ^ M for all positive integral values 
of n. Thus 

( 11 ) 

and also if for all values of | aS' | ^ C, where O' is a constant, 

( 12 ) 


Hence if /(Q)^- ^ a^Q’^S, (13) 

n=0 


each term is less than the corresponding term of the series 

9- 

\r) n\ ’ 


2 M( 

11=0 


.(14) 


which converges for all values of x however large. So long as f{z) is 
expansible about the origin in a convergent power series, however 
small its radius of convergence may be, the expression /(Q) S will be 
the sum of an absolutely convergent series however great x may be. If 
/(;::) and g{z) are both expansible within some circle, their product 
series will also converge within this circle and the expression/(Q) g {Q) S 
will be an absolutely convergent series however great x may be. 
Provided with this result we can now easily extend the result (7) to the 
case where / {z) and g (z) are infinite series, by methods analogous to 
those used to justify the multiplication of two absolutely convergent 
power series. 


1.3. We can now extend these methods to the solution of a family of 
n simultaneous differential e(iuations of the first order with constant 
coefficients. Suppose the equations are 

^ 11^1 +<^ 12 ^ 2 + ••• = ' 

^21 yi + ^22^2 + • •• + e^nyn = N 


^ni^i ^‘n^y^ + • • • + ^nnyn ~~ 
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where the are dependent variables, x the independent variable, era 

denotes ara-^ + Ka. where ar« and are constants, and the are 
ax 

known functions of We do not assume that 

ara — elan ^ra ~ 

but we do assume that the determinant formed by the a’s is not zero. 
When x = 0, yi~Ui, and so on, where the u's are known constants. 
First perform the operation Q on both sides of each equation. We have 


Qer.y, = 




= ar»(ys-U,) + br,Qy, 

= /r«y«-ar.M, (2) 

where /r« denotes the operator + ^i’hen the equations and the 
initial conditions are together equivalent to the equations 

/ii yi +/i2 3^2 + • • • +/in yn^'ei + QSi] 


.(3) 

•( 4 ) 

.(5) 


fni yi +/n2 3/2 + • • • +/nn Vn = t’n + 

where Vr- ariUi+a^u^-^ ... + arnUn. 

The general equation can be written compactly 

^sjra ya ~ Q,^r 

Now let D denote the operational determinant formed by the /'s, 
namely, 

fn fl2 ••• f\n 

f'lx f'n ••• f^n I . . 

•W 


I fnx ful • ’ • fni 
If this determinant is expanded by the ordinary rules of algebra and 
equal powers of Q collected, we shall obtain a polynomial in Q. The 
term independent of Q is simply the determinant formed by the a's, 
which by hypothesis does not vanish. Now let I^ra denote the minor of 
fra in this determinant, taken with its proper sign. jPls also is a poly- 
nomial in Q. 

Now operate on the first of (3) with Fig, on the second with i^ 2 «, and 
so on, and add. Then in the sum the operator acting on ymy say, is 

XFrsfnn (^) 

If m = 8, this sum is the determinant i> ; if m 4= it is a determinant with 
two columns equal and therefore is zero. The resulting equation is 
therefore 


Dya = 2,. Fra (t?r + Q^r) . 


(8) 
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Now if all the are bounded and integrable within the range of values 
of X contemplated, the expression on the right of (8) is also a bounded 
integrable function of Also since the function obtained by 

replacing Q in /) by a number Zy is regular and not zero at z = 0, the 
function l/JD (z) is expressible as a power series in with a finite radius 
of convergence. Define D~^ as the power series in Q obtained by putting Q 
for z in the series for 1/2) (z). Then operate on both sides of (8) with 2)“\ 
We have 

Frs {Vr + Q8r) (9) 

But, since series of positive integral powers of Q can be multiplied 
according to the rules of algebra, D~^D gives simply unity, and we 
have the solution 

ys = D-^XFrs{Vr^QSr) (lO) 


This gives a complete formal solution of the problem. Its form is often 
convenient for actual computation, especially for small values of ; but 
it can also be expressed in finite terms. D is, as we have seen, a poly- 
nomial in Q of degree n at most, while Frs is a polynomial in Q of 
degree w - 1 at most. Our solution is therefore of the form 


„ = tW) + QS 
I)(Q) D(Q)^ 


( 11 ) 


where and xj/ are polynomials whose degree is ordinarily one less than 
that of 2). Since the determinant formed by the a’s is not zero we may 
denote it by Ay and then D will be the product of n linear factors, 
thus 

2) = A (l-'ai0(l-a2^) ... (1 -a„Q), (12) 


where the as will ordinarily be all different. Then ^(Q)/2>(Q) can be 
expressed as the sum of a number of partial fractions of the form 


L 

l-aQ' 


(13) 


But by 1.12 (3) this is the same as The part of the solution arising 
from the -u's can therefore be expressed as a linear combination of 
exponentials. 

The justification of the decomposition of <^(Q)/2)(Q) into partial 
fractions is that this decomposition is a purely algebraic process. Hence 
if the partial fractions and the original operator are all expanded in 
positive powers of Q, and like powers of Q are collected, the coefficients 
of a given power of Q will be the same in both expressions, and the 
equivalence is complete. 
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Exceptional cases will occur if D contains no term in or if two 
or more of the as are equal. If the term in Q" is absent the expansion 
of the operator in partial fractions will usually contain a constant term. 
If the term in is also absent we must divide out, and the expan- 
sion in partial fractions will contain a term in Q; this will give a term 
in X on interpretation. 

If several of the a's are equal, the expression in partial fractions 
will involve terms of the form il/(l These can be interpreted 

by direct expansion ; but another method is more convenient. Starting 
with 

(14) 

let us differentiate r-1 times with regard to a. We find 


( 1 -^ 





(15) 


so that 


Qr-i 

a--Qf 



(16) 


A rather different form of resolution into partial fractions from the 
ordinary one is therefore necessary if each fraction is to give a single 
term in the solution. Instead of having constants in all the numerators 
we must have powers of Q, the power of Q needed in any fraction being 
one less than the degree of the denominator. But if we write p'^ for Q 
the fraction in (16) is algebraically equivalent to pl{p - aj'. In this form 
the power of p in the numerator is independent of the degree of the 
denominator, and it is easier to resolve into partial fractions of this form 
than into those involving Q directly. 

The above remarks apply to the interpretation of the effect of the 
initial conditions; that is, of the first term on the right of (11). To find 
the effect of the terms in the >S’s we may expand the operators simi- 
larly. To interpret (1 - aQ)~^ QS in finite terms, we note that it is the 


solution of ^ = that vanishes with x. This is easily found by the 


ordinary method to be 

y = (17) 


which is the interpretation required. 

To sum up, we can solve a family of equations of the type (1) by 
first integrating each equation once from 0 to ^ with regard to x, allow- 
ing for the initial conditions. This gives a set of equations of the type 
(3). The subsequent process for deducing the operational solution (10) 
is exactly the same as if the operators fra were numbers, and the 
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ordinary rules of algebra were applied. The solution can be evaluated 
by expanding the operators in ascending powers of Q and evaluating 
term by term; this in general gives an infinite series. Alternatively it 
can be obtained by resolving into partial fractions and interpreting each 
fraction separately; this gives an explicit solution in finite terms. 

1 . 4 . Heaviside's method is equivalent to that just given; it differs 
in using another notation, which is not quite so convenient in a formal 
proof of the theorem, but is rather more convenient for actual applica- 
tion. In Heaviside’s notation the operator above called Q is denoted by 
At present we need not specify the meaning of positive powers of 
p ; negative integral powers are defined by induction, so that denotes 
We have noticed that the passage from 1.3 (3) to 1.3 (10) is a 
purely algebraic process. Consequently if all the equations (3) were 
multiplied by constants before the algebraic solution the same answer 
would be obtained. Suppose then that we write the general equation 


1.3 (5) in the form 

^s(j^ra ^rsP )y«“ P ^>S(. . ......(l) 

Multiply throughout by as if this were a constant. We get 

^8 iS^raP pa ~ ^a^raP^s C^) 


On solving the n equations of this form we shall obtain a solution 
identical with 1.3 (10) except that p~^ will appear for Q, and both 
numerator and denominator will be multiplied by the same power of/?. 
If the operators in the solution are expanded in negative powers of />, 
and /?"^ is then interpreted as Q, the result will be identical with that 
already given. Comparing (2) with the original equations 1.3 (1) we see 
that the new form of our rule is as follows : 

Write p for djdx on the left of each equation; to the right of each 
equation add the result of dropping the 6’s on the left and replacing 
the i/’s by their initial values; solve the resulting equations (2) by 
algebra as if p was a number ; and evaluate the result by expanding 
in negative powers of p and interpreting p~^ as the operation of 
integrating from 0 to x. 

This is Heaviside’s rule. In what follows the equations (2) will usually 
be called the subsidiary equations. 

To obtain the solution explicitly, we put era for ciraP + irs, denote our 
determinant 

^11 ^12 ••• 

^21 ^22 • • • ^2?l 


€n\ ^n2 ^nn 
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by A, and denote the minor of in this determinant, taken with its 
proper sign, by Era- Then the solution is 

E 

y* ~ ^ ^r) (3) 


Since the determinant formed by the is not zero, A is of degree n in 
Pj while Era is at most of degree n- 1. The operators can therefore be 
expanded in negative powers of jt?, as we should expect ; positive powers 
do not occur. All terms after the first vanish with iv; the first is 





( 4 ) 


where Ars is the minor of ara in A, But 

XAradrm = A (m = s)\ , . 

= 0 

and (4) reduces to Ua^ as we should expect. This verifies that the solu 
tion satisfies the initial conditions. 


1.5. To interpret the operational solution 1.4 (3) in finite terras we 
require rules for interpreting rational functions of p operating on unity 
and on other functions. We have already had the rules 

= Q ) P~^ = so on, (1) 

p^^l = Ql = x\ p~^l = and in general/?"’*! ...(2) 

If unity is replaced by Heaviside's ‘unit function,’ here denoted by 
H{x\ which is zero for all negative values of x and 1 for all positive 
values, we shall still have 

p~^H{x)^f^, (3) 


when X is positive, but it will vanish when x is negative. We can also 
replace the lower limit of the integrations by - qo without altering this 
interpretation. 

Again, we shall have when x is positive 


P 

p — a 


l—aQ 




(p-af ’ 


( 4 ) 

( 5 ) 


a r) 

= 

p - a p — a 


( 6 ) 


where the function operated on may be either unity or H{x\ In the 
latter case all the operators will give zero for negative values of x. 
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The operators in 1.4 (3) are of the form /(/?)/ where/(jo) and 
F{p) are polynomials in jt?, and f{p) is of the same or lower degree 
than F{p). If F(j>) is of degree n it can be resolved into n linear 
factors of the form p-a. Then provided that the a's are all different 
and none of them zero we have the algebraic identity 

pF(j>) pb\0) "aF' {a) p- a' ^ ^ 

whence ^ + 2. ^ (8) 

F(p) F(0) aF (a)p-a ^ ^ 

If this operates on unity or //(^) we have therefore for positive values 

of ^ 

/Sp) 1 _/(o) ^ ^ ra) 

y{p) FiSty^aF'io.)^ ^ ^ 

To justify this we notice as before that (<S) is a purely algebraic iden- 
tity, and therefore if both sides are expanded in negative powers of p, 
beginning with constant terms, the expansions of the two sides will be 
identical, and on interpretation in terms of integrations will give the 
same result. 

The formula (9) is usually known as Heaviside’s expansion theorem ; 
but as Heaviside’s methods involve two other expansion theorems* it 
will be called the ' partial-fraction rule ’ in the present work. 

If some of the a’s are equal or zero, the expression (7) considered as 
a function of p will have a multiple pole, and its expression in partial 
fractions will contain terras of the form (;> - a)"*, where s is an integer 
greater than unity, and a maybe zero. Then will contain 
terms of the forms or pj^p — a)% which can be interpreted by 

means of (2) or (5). 

By means of these rules we can evaluate all the expressions for the 
part of in 1.4 (3) that depends on the initial values of the yB. If the 
^’s are constants for positive values of x, as they often are, the same 
rules will apply to the part of the solution depending on them. If they 
are exponential functions such as the easiest plan is usually to 
rewrite this asp/(/?-/A) and reinterpret. Thus 

F{p)^' F{p)p-y- ^ 

= ^ + 2 y — y eT ( 11 ) 

F(,x) (a-M)i' (a) 


* Namely, expansion in powers of Q or p~^ and interpretation term by term; 
and expansion in powers of e~P\ where /i is a constant, as in 4.2. 
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If S is expressed as a linear combination of exponentials we can apply 
this rule to each separately. This is applicable to practically all func- 
tions known to physics. 

Alternatively we can resolve the operator acting on S into partial 
fractions and interpret p~'^S by integration and other fractions by the 
rule of 1.3 (17) 

-L-8 = ^ rSe-'^dv ( 12 ) 

p-a Jo 

This completes our rules for solving a set of linear equations of the 
first order with constant coefficients. In comparison with the ordinary 
method, we notice that the rules are direct and lead immediately to a 
solution involving operators, which can then be evaluated completely 
by known rules. If it happens that we only require the variation of one 
unknown explicitly, we need not interpret the solutions for the others. 
In the ordinary method we have to find a complementary function and 
a particular integral separately. To find the former we assume a solu- 
tion of the form and on substituting in tlie differential equa- 

tions with the omitted we find an equation of consistency to 
determine the n possible values of a, and the ratios of the X’s corre- 
sponding to each. The particular integral is then found by some method, 
but it does not as a rule vanish with x. The actual values of the Vs 
are still undetermined, and tlie value associated with each a must be 
found by substituting in the initial conditions and again solving a set 
of n simultaneous equations. The labour of finding the equation for the 
a’s and the ratios of the A's corresponding to one of them is about the 
same as that of finding the operational solution in Heaviside's method ; 
the rest of the work is avoided by the operational method. Further, if 
some of the a's are equal or zero considerable complications are intro- 
duced into the ordinary method, but not into the operational one. 

1.6. The above work is applicable to all cases where A, the deter- 
minant formed by the a's, is not zero. We can show that if A is zero 
there is some defect in the specification of the system. A system is 
adequately specified if when we know the values of the dependent 
variables and the external disturbances the rates of change of the 
dependent variables are all determinate, and if the initial values of the 
variables are independent of one another. Now if A is zero, let us 
multiply the r’th e(iuation by Ars, the minor of in A, and add up 
for all values of r. The coefficient of dygjdx in the sum is A, which is 
zero ; that of any other derivative is a determinant with two columns 
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equal, and is also zero. Thus the derivatives disappear entirely from 
the sum, and we are left with a relation between the ^s and the 
If the coefficients of the also vanish, there must be a permanent 
relation between the S's, and one of the equations we started with is a 
mere logical consequence of the others; then we have not enough 
equations to determine the derivatives. If the coefficients of the ^s do 
not vanish, we can put and obtain a relation between the initial 
values, which are therefore not independent. 


1.7. The method is most easily extended to equations of higher 
order by breaking them up into equations of the first order. Thus if we 
have an equation of the second order such as 


dhi dy , 


( 1 ) 


we introduce a new variable given by 



( 2 ) 


and the original equation can be replaced by 


dz 

= aS’. 


( 3 ) 


We have now two equations of the first order in y and z. If initially 
y=y<, and dyjdx-yi, the subsidiary equations are 


py-z=pyo, (4) 

{p + a)z + by=p^ji + S. (5) 

Solving by algebra we find 


(/ + ap) //o + pi p + < S* 
p- + ap+b 


To interpret, put 

p‘ + np + h = {p-o){p-p), 
and apply tlie partial-fraction rule. We find 


( 6 ) 

( 7 ) 


y = ^ ^ ^ (2/i - ^y.) ^ 


+ j’^'se-P-^dx-e^^IJse-'^dxJ, ...( 8 ) 

which is easily shown to satisfy all the conditions and therefore to be 
the solution required. 
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1 . 71 . A few illustrative examples may be given, {p is written for 
djdx from the start.) 

1. + + 3)y = l; yo=3; yi = -2. 

We consider the subsidiary equation 

+ 3) y = 1 - + 3 (/)“ + ip) 

= 3/ + 10y+l, 

Zf + l Ojp + 1 
^~\p+l){p + Z) 

= J + 3e-*-Je-”, 

on interpreting by the partial-fraction rule. 


2 . (p’-i-5j» + 6)y = 12; yo=2; y, = 0. 

Consider (p’ h- 5jo -i- 6) y = 12 + 2 (p’ + bp). 

Cancelling the common factor, 

y = 2. 

3. (p + 3)y = e-’^; y(, = 0. 

This can be written 

_ p 
y (p + 2)(p + 3) 

{p^ + y = 0. 

y<s + Vy\ sin nx 

y ~ yo cosnx + pi 




or 


4 . 


6. (p + 2)=y = a^e-“®; y„ = 0; y, = 0. 

The expression on the right is equivalent to the operator 2p/(p -t- 2)’. 


Hence 


2p 2ar^ 


12 


x^e~ 


We notice the advantage of Heaviside^s method in avoiding the use 
of simultaneous equations to determine the so-called ‘arbitrary ^ constants 
of the ordinary method. In particular in the second example the data 
have a property leading to a simple solution. Heaviside's method seizes 
upon this immediately and gives the solution in one line; with the 
ordinary method the simultaneous equations for the constants would 
have to be solved as usual. 
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In general we may say that the more specific the problem the greater 
will be the convenience of the operational method. 


1 . 72 . The method just given can be extended easily to a family of 
equations of the second order. If the typical equation is again written 



2^ Crg yg — Sr , 

(1) 

where now 

, d 

er.-ar.^^2+<>r,^^ + C„, 

(2) 

we can introduce 

n new variables Zi^ given by 



1 

II 

O 

(3) 


thus treating the first derivatives of the y’s as a set of new variables. 
Then the equations (1) are equivalent to 


(^r« ^8 + y«) ~ (4) 

and (3) and (4) constitute a set of 2;^ equations of the first order. Suppose 
also that when ^ = 0, = According to our rule of 1.5 we 

must replace djdx by p and add pUr to the right of (3), and 2, arspvg 
to the right of (4). We have now to solve by algebra, and may begin by 
writing the revised form of (3) 

Zr^p{yr-Ur) (5) 

When we substitute in the modified form of (4) we get 

X (lr»P^ (j/s - Us) + 2 , brsP (Ps - Ug) + 2 , Crg + 2 * arg pVg, ( 6 ) 

or, on rearranging, 

2. (arsP^ + brsp + Crg) ^g = 2, P^ + p) Ug + 2« Urg pV, + . . . . (7) 

We have thus n equations for the y's to solve by algebra, and the 
solution can be interpreted by the rules. 

In Bromwich's paper ‘Normal Coordinates in Dynamical Systems' a 
method equivalent to the operational one is applied directly to a set of 
second order equations. First order equations, however, arise in some 
problems of physical interest, and merit a direct discussion. This, as we 
have seen, is easily generalized to equations of the second order. 

1 . 73 . In the discussion of the oscillations of stable dynamical systems 
by the ordinary method of normal coordinates there is a difficulty when 
the determinantal equation for the periods has equal roots. This does 
not arise in the present method. If the system is not dissipative, 
and the determinant formed by the ^'s as defined in (2) has a multiple 
factor {p^ + we know from the theory of determinants that every 
first minor has a factor (p^ + On evaluating the operational solution 
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of (7) we can therefore cancel the factor from the numerator 

and the denominator, and we are left with a single factor + a^) in tlie 
denominator. Thus the part of the solution depending on the initial 
conditions is of the trigonometric form as usual ; terms of the forms 
t cos at and t sin at do not arise. 


1 . 8 . In all the problems considered so far the operators that occur in 
the solutions are expansible in positive powers of Q, or in negative powers 
of p, and we have seen that so long as the series obtained by replacing 
Q by a number have a finite radius of convergence, however small, the 
operational solution is intelligible in terms of the definitions we have 
had. But we have not defined p as such, because we have only needed 
to define its negative powers, and p cannot be expanded in terms of its 
own negative powers. Then has p any meaning of its own? Since it 
replaces djdx when we form the subsidiary equation we may naturally 
suppose that/) means djdx, and this is the meaning sometimes actually 
attributed to it; but care is needed. We recall that when the subsidiary 
equation is formed a term like py^ appears on the right ; but dy^jdx is 
zero, so that if we pushed this interpretation too far we should be faced 
with the alarming result that the solutions of the e(iuations do not 
depend on their initial values. The fact is that though the operators 
djdx and Q both satisfy the laws of algebra and are freely commutative 
with constants, they are not as a rule commutative with each other. 
Thus j X rx 




^ 

Thus the operators p and Q are commutative if, and only if, the function 
operated on vanishes with x. It appears from (1) that d\dx undoes the 
operation Q if dldx acts after With this convention we can identify 
/), the inverse of ft with djdx. When p and Q both occur in an operator, 
the Q operations must be carried out before the differentiations*. 

This result explains why some of our interpretations differ from those 
given in text-books of differential equations for the determination of 
particular integrals. For instance, we have interpreted l/(/)~a) as 

- {e^ - 1). In the ordinary method, due to Boole, we expand this operator 
in ascending powers of /), and get 


1 

p- a 


1 - 




* Cf. Heaviside, Electromagnetic Theory^ 2, 298. 


( 3 ) 
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This assumes that the sum of all the terms of the series after the first 
is zero. But 


( 1 -. 1 ) 1 = 1. 

\p- a a/ a [p — a) 


.(4) 


If the operator on the right is written clearly 


zero. 


But if it is written - . jt? . . 1, it gives by our rules - The differ- 

a ^ p~ a ^ a 

ence arises from the fact that the operators p and {p — a)"^ are non- 

commutative. 

We may notice, incidentally, that whereas the series in powers of Q 
give convergent series on interpretation, the same is not true of the 
corresponding series in powers of p. For instance, \i S-{x- I)'", where r 


is fractional, the series for 

p-\ 


S in ascending powers of p diverges 


like 2 Apart from the greater internal consistency of 

Heaviside’s methods, then, they are capable of much wider application. 
We have also 


eVC-r) = + ^ + • • •)/(a;) 


=/W + A/'(^) + |V"W+- 


=/(^ + A), (5) 

by Taylor’s theorem. The operator increases the argument of a 
function by h. 

This operator is useful in enabling us to find a theorem that plays 
in Heaviside’s method a part closely analogous to that played by Fourier’s 
integral theorem in the ordinary method. Instead of trigonometric 
functions, Heaviside treats as fundamental the function here called 
II{x)y which is zero for negative, and unity for positive, values of 
This function being discontinuous, we cannot at once apply Taylor’s 
theorem to it; but continuous functions can be found as near to it as 
we liket, and the results obtiiined by applying (5) may be regarded as 
the limits of results proved for these. Then we shall write 

= 1 when x >- /n 
0 when x <- i 


J 


* Proc, Poy, Soc,^ A, 52, 513. 

f E.g. 4(Erf>u;-f 1) where X may be made indefinitely great. 


2 
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A fuller justification of this step will be given in the next chapter. Then 
if hi < hf2 > 

H{x)-\ if hi<x< Jh 1 

=:0if;r<Ai i (7) 

and = 0 if .^* > ^2 J 

If further Aj < As < ^3 < . . . < 

+ • • • /(An)} /^(^) 

= 0 for ^ < Ai , /(hi) for hi <x < A., and so on (8) 

If then the subdivisions of the interval fh to An become indefinitely 
numerous, we can approximate as closely as we like to a function f (x) 
by a sum 2 / (A) H (x). Now with a further proviso about 

replacing H(x) by a continuous function and later proceeding to a 
limit, we can replace this sum by an integral. Then 

/(/) = r f(h) d H(x)] - r f(h) pH(x) dh. (9) 

J h = —OO J -CO 

If the integration with regard to A is carried out, we obtain a function 
of p operating on H (x). In this way a given function of x can be 
expressed in the operational form*. 

We notice that if A is positive 


e-P^- H(x) = eri>’^ 
V 


( 0 when x < 0 
1 X when x > 0 


{ 


0 when x<h 
x-h when x > 


h, 


and 


^ e-p>‘ 11 {x) H (a: - A) = 11 {x - h) dx 

~ 0 when x<h. 


-x-h when x > A. 

Thus the operators ^ and can be commuted provided A is positive. 

1 rr / X f 0 when x <-h 
But ev^-H(x)^{ 

p ^ \X‘ 


and 


i.f- //«={" 


+ A when x>-h, 
0 when x < 0 
; when X > 0. 


Fortunately symbolic solutions seldom contain operators of the form 
where A is positive, so that the fact that this operator is not com- 
mutative with - gives little trouble in practice. 


Heaviside, Electromagnetic Theory^ 3, 327. 
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CHAPTER II 

COMPLEX THEORY 

2.1. All the interpretations found so far for the results of operations 
on unity or on II {ar) are special cases of two closely related general 
rules, as follows. If </> is an analytic function, 

( 1 ) 

= ( 2 ) 

In the former definition the integration extends around a large circle 
in the complex plane. In the latter the integration is along a curve 
from c — too toc + tcc, where c is positive and finite, such that all 
singularities of the integrand are on the left side of the path (that is, 
the side including - x). These rules are botli due to Bromwich*. 
Considering first (1), we have 

■;■■■ w 

and the only pole of the integrand is at the origin, where the residue 
is !. Hence when 7i is a positive integer 

.,n 

« 

Also we see that 

(5) 

since the integrand has no singularity in the finite part of the plane ; 
and 

; ( 6 ) 

The interpretation 1.5 (5) of also follows immediately from 

The partial-fraction rule can also be proved easily. For by (1), with 
the notation of 1.5, 

/M ZW (7) 

_ , / ( 0 ) / ('*) ax 

«>’'(«) 

on evaluating the residues at the poles 0, ag, 

^ ‘Normal Coordinates in Dynamical Systems,’ Proc. Lond. Math. Soc., 15 (1916), 
401-448. 


2-2 
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We can also prove that 1.4 (3), when interpreted by the rule (1), is 
the correct solution of eiiuations 1.3 (1). If we denote by e„(K), A (k) 
and £'r, (k) the results of replacing jp by « in e,.», A, and the part 
of 1.4 (3) not involving the (S’s is equivalent to 






(k) __ J 

A (,) 


.(9) 


Substituting in the differential equations 1.3 (1) we find that the left 
side of the w’th equation is 


1 


But 


a"m (10) 

A (k) if r = m 

(11) 


27rc./e^"~“ A(<c) 

2s i/rs ('<) (x) = A («) if r = m 

= 0 if r 4= 7n, 

and the expression (10) is then equal to 




.( 12 ) 


Thus the differential equations are satisfied. 

Now considering the initial conditions, we put and find 



1 /* ^ PrsCO 7 

am”'* 

(13) 

Also 

'^r ^ {^nn (^) ^rm\ ^ud 

(14) 

and 

/c K A (k) ““ 

" ^ K K a’'(k) 

(15) 


Now A (k) is a polynomial of the w'th degree in k, and (k) one 
of the (n - l)th degree. When | k j is great enough the second term in the 
integrand is of order It therefore gives zero on integration, and we 
have 

(Iff) 

so that the initial conditions are satisfied. 

Next, suppose that instead of the>S"s being zero they are exponential 
terms of the form Since the solutions are additive it is enough 

to suppose that the initial values of the j/’h are all zero, and to consider 
the effect of an exponential term in the first equation alone. Then the 
equations are 

^rs 2^8 “ 1 ) 


••( 17 ) 



COMPLEX THEORY 


21 


The operational solution is 


This is to be interpreted as 


y» = 


Eu Pp 

^ p-li' 


y> 


1 f Eu (k) 

2vi}c A ((f) 


d/c. 

K-fX 


(18) 

(19) 


Substituting in the differential equations we find that the left sides 
of all vanish except the first, in consequence of the integrand con- 
taining as a factor a detenninant with two rows identical; the first 
gives 


Zttl Jfj k — 


( 20 ) 


The differential equations are therefore satisfied. 

When and \k \ is large, the integrands in (19) are all 0 (k~^) and 
the integrals are therefore zero. Thus all the vanish with a:. 

Since all functions that occur in physics can be expressed by Fourier s 
theorems as linear combinations of exponentials, the proof that the 
solution can be carried out by operational methods is complete. We 
shall see later, however, that Fourier's theorems are not so convenient 
to use as the formula 1.8 (9) based on the function //(.r). 


2 . 2 . Now let us consider the integrals 2. 1 (2), taken along the path L, 
If we suppose the contour completed by a large semicircle from c + too 
to c- too by way of — x , all the singularities of the integrand are within 
the contour, and therefore tlie integral around it is the same as the 
integral (1) around a large circle. Now if x is positive, and {k)/k tends 
to zero uniformly with regard to arg k as k tends to x , the integral 
around a large semicircle tends to zero as the radius becomes indefin- 
itely large, by Jordan’s lemma^. Hence the integral alongZ is equivalent 
to the integral around a large circle. Thus if x is positive and 

(ac)/k - 0 (k-’») 

when K is great and n is positive, 7/ (.r) is the same as 1. 
But if X is negative, the integral around the large semicircle on the 
negative side of the imaginary axis is no longer an instance of Jordan's 
lemma, and this result no longer holds. The integral around a large 
semicircle on the positive side of the imaginary axis is, however, re- 
ducible to the integral along Z, since by construction there is no 
singularity between these two paths. If then <^(k)/k = 0 (k"”) when k is 


* Whittaker and Watson, Modern Analyeiz^ 1915, 116. 
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great, the limit of the integral around this large semicircle is zero by 
Jordan's lemma, and therefore ^ (p) II (^r) is zero. 

Thus if (p) is expansible in descending powers of p, beginning with 
a constant or a negative power of p, we shall have 

(p) H {x) = <#> (jt?) 1 when > O') 

<l){p) II {x)-0 when.^<0/ 

The integrals 2.1 (2) for p'^H{x\ where ?? is a positive integer, are 
divergent, but these derivatives can be found by differentiation. Evi- 
dently they are zero except when .r = 0, when they do not exist. 

Again, (u') = dK = H{x + h) (2) 

This proves the result obtained by less satisfactory means in 1.8. 

We can also translate into the form of a double integral the opera- 
tional expression for a general function. From 1.8 (9) 

f{x) = f (3) 

J -‘00 

where the integration* gives a function of ^ operating on II (x). Inter- 
preting by Bromwich’s rule 

/W = l/W (4) 

where the integration with regard to ^ is to be carried out first. 

If the function of p that arises when (3) is integrated has no singu- 
larity on the imaginary axis or on the positive side of it, and if further it 
contains p as a factor, we can replace the integration with regard to k 
in (4) by one along the imaginary axis. Then we can put 


K=:lX, ( 5 ) 

and we have 

/(x)~~ I f /(A) sinX(x-/^) MdX (6) 

TT J-CC Jo 

This is Fourier’s integral theorem. 


2.3. In practice the form 2.1 (2) is generally used in preference to 
2.1 (1). In problems involving a finite number of ordinary differential 
equations of the first order the forms are equivalent for all positive 
values of the independent variable ; and since the independent variable 
is usually the time, and we require to know how a system will behave if 
it is in a given state at time 0 and is afterwards acted upon by known 

* The integral is of the type introduced by Stieltjes. 
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disturbances, it is as a rule only positive values of the time that con- 
cern us. But when we come to deal with continuous systems it usually 
happens that when the operational solution (f> {p) is interpreted as an 
integral, the integrand has an infinite number of poles, and that no 
circle however large can include all of them. Hence the contour integral 
2.1 (1) cannot be formed. But the line integral 2.1 (2) usually still 
exists. Again, it may happen that (k) has a branch-point at the origin 
or on the negative side of the imaginary axis. Here again the contour 
integral does not exist, but the line integral does. 

Conse(piently in most of the writings of Heaviside and Bromwich, 
when a function of an operator (/?) occurs without the operand being 
stated explicitly, it is to be understood that the operand is H{x\ and 
that the interpretation 2.1 (2) is to be adopted. This rule will be 
followed in this work when we come to treat continuous systems; and 
if it is also supposed to be adopted in the problems of the next chapter 
no harm will be done. 

Considerable latitude is admissible in choosing the path L. Since by 
construction the integrand is regular at all points to the right of Z, we 
may replace L by any other path on its positive side provided it ends 
at c ± t X . Again, the systems we are to discuss are usually stable 
systems ; that is, the poles of <\> (k) are all on the imaginary axis or to 
the left of it. Then Z can be taken as near as we like to the imaginary 
axis. Z must not cross the imaginary axis if there are poles on the 
latter; but it can be taken to be a line parallel to the imaginary axis 
and distant c from it. This is the device most often used by Bromwich; 
the only advantage of the present definition of Z is that the solution 
is then adaptable to unstable systems as well as stable ones. In all cases 
the actual value of c does not aflect the results, so long as it is 
positive. 


2 . 4 . The principal operators involving branch-points are fractional 
powers of p. We re(iuire then an interpretation of p^^H (.r), where n is 
fractional*. By our rule 

//*Zr (.r) = ^ I (1) 

and the integral converges if n < 1. A contour including Z as part of 
itself, and such that the integrand is regular within it, is as shown. 
Evidently if x is positive the large quadrants make no contribution. 


Cf. Bromwich, Proc, Camh. Phil. Soc.y 20 (1921), 411-427. 
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The integral along L{= A B) is therefore in the limit equal and opposite 
to one along CDEF, If further n is positive the contribution from DE 
tends to zero with the radius of DE, Thus A B, CD and contribute 
the whole of the integral. Now on Cl) and EF 

K~ and /x (2) 

respectively, where /x is real and positive. Hence on CD 

and on \ (3) 

= - /x"-^ erf^^ dfi ) 



Fig. 1. 


Therefore 

1 /■«> 

p^H {x) ~ — I /Ex”“' dfx sin (w - 1) TT 
_ sin mr r {n^ 

TT 


But by a known theorem 

r (t^) r (1 ~ ^^) = TT cosec WTT, 


(4) 

(5) 


fH{x) = 


x-^ 


and therefore 


( 6 ) 
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This result has been obtained for {}<n<l. If we change n into -tw, 
where “ 1 < w < 0, we have 


^ ^ ^ r(w + 1) 


Powers of p outside this range can be found now by integration or 
differentiation. We have 

{x) > (8) 


r + \\my 

which shows that (7) still holds when m- \ is written for m. (To justify 
this, tlie path AB must be replaced by FEDC.) 

Similarly 

p-^^n{x) = tn ~ ( 9 ) 

yol O^^^ + l) r(;« + 2) ^ ^ 

By induction we may therefore generalize (7) to any value of m. Hence 

O") 

Since when m is a positive integer r(m + l) = 7?i !, and when m is a 
negative integer V (771 + 1) is infinite, the interpretations already adopted 
for integral powers of p are special cases of these. 

In particular, since 

r(|)=V^, (11) 

we have 

X 1 5 1 -2 ^ 1 I -5 

“i (1^) 


and so on. Also 


) 




A related function that arises in problems of heat conduction is 
where a is a constant, and q denotes p^. By Bromwich’s rule 

1 r pKX - aA 

. -* ('« 

On L the argument of is between ± Hence if a is positive the 
integral is convergent. 

Immediate expansion of in a power series and integration term 
by term would not be legitimate, because all the resulting integrals 
after the first two would diverge. But (14) is equivalent to an integral 
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along a path such as FEDC in Fig. 1, and on this path we can proceed 
in this way. Thus 




1 a^K 

aK^ +21“ 


...(15) 


All the positive integral powers of k give zero on integration. The other 
terras are equivalent to 

2 t » _ 1/ « ^ 

3 ! 5 ! n/tt l2 J.V 3 V2 2 \ 5\2 Jx / J 


1-Erf 


where, by definition, 


- / (17) 


By differentiation with regard to a we find 

( 18 ) 

We shall sometimes need an asymptotic approximation to Erf w when 
w is great. We have 

1 - Erf w = I f ir^ da 

TT J tv \> TT J 

1 , , r. 1 . 1.8 , 1 . 3.5 , "1 

= IV 1 - - v'~~ + — iv~* — (t’~* + . . . 

Vtt L 2 - 2.2 - 2 . 2.2 J 


( ll I rJ ) 

2. 2. 2... 2 


C”’* ir^ ( 2 n+l) (19) 


on successive integrations by parts. But the last integral 

=-2 [ d^ < f 


< ar (2« - 0 (20) 

2n-l ^ ^ 

SO that its contribution to the function is less than the previous term. 
We have therefore the asymptotic approximation 


^ (1 - Erf w) ~ 

s/tt L 2 


3 1.3 3 1.3.5 

^,-J ^ ^ 3 w ^ ^ 


2 2.2 2 . 2.2 


..]. ( 21 ) 
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CHAPTER III 


PHYSICAL APPLICATIONS: ONE INDEPENDENT VARIABLE 


3 . 1 . An electric circuit contains a cell, a condenser, and a coil with 
self-induction and resistance. Initially the circuit is open. It is suddenly 
completed ; find how the charge on the plates varies with the time. 

Let y be the charge on the plates, t the time, K the capacity of the 
condenser, L the self-induction, and R the resistance of the circuit; 
and let E be the electromotive force of the cell. Write or for djdt. 

The current in the circuit is y, and the charging of the plates of the 
condenser produces a potential difference yjK tending to oppose the 
original E.m.f. Then y satisfies the ditferential equation 

+ ( 1 ) 


Initially y and y, the 
is simply 


and the solution is 


current, are zero. Hence the subsidiary equation 


\ A 


E 

y~^ 

Zcr^ -I- R(t 



If now the denominator is expressed in the form A ((t -e a) (o- -i- ^), the 
interpretation is, by 1.5 (9) 



KL + 


1 

a-P 



(4) 


Since a + p and are both positive, a and /? must either be both real 
and positive, or conjugate imaginaries witli positive real parts. In either 
case y tends to KE as a limit, as we should expect. 

We notice incidentally that if the circuit contained no capacity or 
self-induction the differential equation would be simply 


Rny^E, (5) 

Hence if a problem has been solved for simple resistances, self-induction 

and capacity can be allowed for by writing Act + A + for R, For this 

reason this expression is sometimes called a ‘resistance operator,’ and 
the method generally the ‘method of resistance operators.’ 
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3,2. The Wheatstone Bridge method of determining Self-Inductim. 
In this method the unknown inductance is placed in the first arm of a 
Wheatstone bridge ; the fourth arm is shunted, a known capacity being 
placed in the shunt*. 

First consider the ordinary Wheatstone bridge, the resistances of the 
arms being Ri, R^, B^, let x be the current in B^ , y that vt\R^,g that 
through the galvanometer, and h the resistance of the battery and leads. 



Then R^x- R^y + Gg = 0 , 

RiX - Riy -{R^^ Ri^ G)g = 0, 
h{x + y)-\- IL.y + Ri{y + g) = E, 
and on solving we find 


( 1 ) 

( 2 ) 

( 3 ) 


9 


x + y 


Bi R3 ~ Ri Ri G (/fi + Ri + R3 + Rf) + (/ij + B.^ (^Rs + R^ 

If b is large compared with R^ and B 4 , we have nearly 

x + y^Ejb, 

and « = { Elb){R.Ji3-B,R4) 

^ G{B4 + R, + R3 + B4) + {R, + /y {B3 + i? 4 ) ■ ■ ■ ■ ■ 

According to the result of 3.1, we can allow for the self-induction in 
the first arm by replacing i?, by Lx + R^. Let the arrangement in the 
fourth arm be as shown (Pig. 3). The resistance of the main wire is R^, 


( 4 ) 

.(5) 

.( 6 ) 



Pig. 3. 

* Bromwich, Fhil. Mag. 37 (1919), 407^419. Further references are given. 
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that of the shunted portion of it r. Suppose the shunt to have a 
resistance S. Then the effective resistance of the whole arm would be 


P rS „ 

r + S' 


( 7 ) 


If the shunt contains a capacity K, we must replace ^Sby S-^l/Ka-. 
Hence in the formula for g we must replace Ri by Lcr + and R 4 by 


R4- 


r^Kcr 


(r + S) K(t 4- 1 ’ 


The result expresses the current through the galvanometer when the 
battery circuit is suddenly closed. 

It can be shown that in actual conditions g cannot vanish for all 
values of the time. A sufficient condition for this would be that the 
operator R^R^ - R 1 R 4 should be identically zero; then g would be 
identically zero whatever the remaining factor might represent. But 
with our modifications this factor becomes 


( 8 ) 

Multiplying up and equating coefficients of powers of cr to zero, we 
find 


R 4 (r^S)^r\ (9) 

- LR 4 + (/4 R, - Ri R,Xr + S) K+ R, 0, (10) 

74/4 -^1^4-0 (11) 


From the construction of the apparatus r ^ R^, r ^ r + S\ Thus (9) can 
hold only if r = JK4 and S = 0 , The shunt wire must be attached to the 
ends of R 4 , and must have zero resistance. Substituting in (10) we find 


L^R^R^K, ( 12 ) 

together with the usual condition (11) for permanent balance of the 
bridge. 

Actually these conditions for complete balance cannot be completely 
satisfied ; but for the determination of L it is not necessary that they 
should. Suppose the galvanometer is a ballistic one, and that it is so 
adjusted that there are no permanent current and no throw on closing 
the circuit. Thus 


Lim^=0; gdt^^O (13) 

^-^ao J 0 


If g is expressed in the form/((r)/i^^(cr), we have 
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where all the a^s, in the conditions of the problem, will have negative 
real parts. Then the condition that g shall tend to zero gives 


/(0)=0 (15) 

Also j = ^ (16) 


Equation (15) shows that the operational form of g contains o- as a 
factor. Also we can write 




./W 


.(17) 


afi'* (a) cr — a 

SO that (16) is the limit of (- ^/o-) when o- tends to zero. The vanishing 

of Lim g and j gdt imply therefore that the operational form of g 

contains o-- as a factor. Hence the modified form of - RiRa, must 
contain as a factor. Then (10) and (11) still hold; (9) no longer 
holds. We now have 


RJ},-^R,R,=^0, (18) 

LR^~ RiT^K, (19) 


which gives the required rule for finding L. 


3.3. The Beumograph. In principle most seismographs are Euler 
pendulums — pendulums with supports rigidly attached to the earth, so 
that when the earth’s surface moves it displaces the point of support 
horizontally and disturbs the pendulum. The seismograph differs from 
the Euler pendulum as considered in text-books of dynamics in two 
ways : instead of being free to vibrate in a vertical plane, it is con- 
strained to swing, like a gate, about an axis nearly, but not quite, 
vertical, so that the period is much lengthened ; and fluid viscosity or 
electromagnetic damping is introduced to give a frictional term pro- 
portional to the velocity. The displacement of the mass with regard to 
the earth then satisfies an equation of the form 

x-k- 2 kx + n^x = \^ ( 1 ) 

where ^ is the displacement of the ground, and k, n, and A are constants 
of the instrument*. Suppose first that the ground suddenly acquires 
a finite velocity, say unity. Then 

( 2 ) 

* Some instruments, such as that of Wiechert, are not on the principle of the 
Euler pendulum, but nevertheless give an equation of this form. 
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and our 

subsidiary equation is 



(cr^ + 2k(T + n^)x= korll (t) 

(3) 

Put 

cr^ 4- 2#c(r + = (or + a) (a + (3) 

(4) 

Then 

X-, ..Hit) 

(<r 4* a) (o* + /3) ^ 

(5) 


= 0 when ^<0 

(6) 

and 

“ — ^ when ^>0 

CL p 

(7) 


The recorded displacement x therefore begins by increasing at a finite 

1 


rate A, reaches a maximum X 




-18 


after a time 


a-^ 


and 


then tends asymptotically to zero. 

If a and /5 are real, and 13 less than a, we see that the behaviour after 
a long time depends mainly on e~^^; now as the experimental ideal is 
to confine the effects of a disturbance to as short an interval afterwards 
as possible, we see that we should make p as large as possible. But 


p-x-ix'-n-y— ;; - , 

/C4-(K--n^)- 

(8) 

and for a given n, (3 is greatest when k ^ n. This is the condition for 
what is called aperiodicity ; the roots of the period equation are equal, 

real, and negative. Many seismographs are arranged so as 
condition. The solution is then 

to satisfy this 

x=.^^Hit) 

(o- + nj ^ ^ 

(9) 

= 0wlien^<() 

(10) 

-Xte~^^ when ^>0 

(11) 

The maximum displacement is now at time 1/n after the start, and is 

equal to X/en- 


If K<7i, we can put 


1 

II 

(12) 

Then the solution is 


,r - - e”** sin yt, 

y 

(13) 
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and the motion does not die down so rapidly as for K~n, The aperiodic 
state therefore gives the least motion after a long interval for a given 
value of n. 

In practice, however, k is usually made rather less than n; in the 
Milne-Shaw machine, for instance, k is about O'ln. The motion is then 
oscillatory, but the ratio of the first swing to the second is or about 
20. But a* vanishes after an interval Trjy from the start, or about 4//^, and 
ever afterwards is a small fraction of its first maximum. The reduced 
damping effect after a long time is considered less important than the 
quick and complete recovery after the first maximum. The time of the 
first maximum is 1*1/?^ from the start, as against 1/71 for the aperiodic 
instrument and 1*57/;^ for the undamped one. 



3.31. The Galitzin seismograph is similarly arranged, but the motion 
of the pendulum generates by electromagnetic induction a current, 
which passes through a galvanometer. If x is the displacement of the 
pendulum, and y that of the galvanometer, the differential equations 
are 

+ ( 1 ) 

y 2 k, (2) 

Supposing both pendulum and galvanometer to start from rest, we have 



( o ^ f 2k^(t + Ui) (&^ + 2k^(t + n^) ’ 


( 3 ) 
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As a rule the two interacting systems are so arranged that k and n are 
the same for both, and both are aperiodic, so that <c = w. Then 


y= 


(o- + nf 


(4) 



If we suppose as before that 

i = //(<). (5) 

we have y = y H {t), 

= Xm ™ (1^ er"') when t > 0, 

= ( 6 ) 


The indicator therefore begins to move with a finite acceleration, instead 
of with a finite velocity as for the pendulum. The maximum displace- 
ment follows after time (3 - sjs)ln- 1*27/;^, the mirror passes through 
the equilibrium position after time 3/w, and there is a maximum dis- 
placement in the opposite direction after time 4*73/7^. The mirror then 
returns asymptotically to the position of equilibrium. The ratio of the 

two maxima is + \/3)^ = 2’3. In comparison with an instrument 

such as the Milne-Shaw, recording the displacement of the pendulum 
directly, the Galitzin machine gives the first maximum a little later, 
the first zero a little earlier, and the next maximum displacement is 
larger in comparison with the first maximum. 


J 


3 
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In an actual earthquake the velocity of the ground is annulled by 
other waves arriving later ; the complete motion of the seismograph is 
a combination of those given by the separate displacements of the ground. 


3.4. Resonance. A simple pendulum, origiiiallyhanging in equilibrium, 
is disturbed for a finite time by a force varying harmonically in a period 
equal to the free period of the pendulum. Find the motion after the 
force is removed. 

The differential equation is 


x + n^x -/sin nt, (1) 

( 2 ) 

(T- + ^ ^ 

. 7i(T 

The solution is then ^=/ 7 ^ (3) 

+ n^y ^ ^ 

nothing having to be added on account of the initial conditions. 

To evaluate this, we notice that 


o* 


- sin 7it. 
n 


Differentiating with regard to n, we have 


and 


2ncr 

(cr^ + n^y 


, = - ^ cos nt — „ sin nt. 

' 71 

(sin nt - nt cos nt) 


(4) 

(5) 

( 6 ) 


Suppose the disturbance acts for time r7^//^, where r is an integer. At 
the end of this time 




(7) 


The subsequent motion is therefore given by 

= -^cos«!! (8) 


3.6. Three particles of masses w, and tw, in order, are attached 
to a light stretched string, the ends of the string being fixed. One of the 
particles of mass m is struck by a transverse impulse /. Find the sub- 
sequent motion of the middle particle. (Intercollegiate Examination, 1923.) 

If Xi, X 2 y x-i are the displacements of the three particles, P the tension, 
and I the distance between consecutive particles, we find in the usual 
way the three equations of motion 



ONE INDEPENDENT VARIABLE 


35 



I 

1 

II 

( 1 ) 


I = - A (- iTi + 2^2 - iTs), 

( 2 ) 


X3 = ~\{-Xi + 2x^, 

(3) 

where 

11 

(4) 

Initially all the displacements are zero, ^ 2=^3 = 0 , = Ijm. Hence the 

subsidiary equations are 


( 0 ^ + 2 X) Xi - A ^2 = O-Z/w, 

(5) 


- A^i + (l^cr^ + 2 X) X 2 - Xxs = 0 , 

( 6 ) 


- \X2 + (or^ + 2 X) Xs = 0 

(7) 

Hence 


( 8 ) 

and 

/21 „ . 2\^ \ Act / 

\2() 0-2 + 2A/ ^ (7^ 4 - 2A m 

(9) 

Multiplying by o-^ + 2A. we find 


and 

(70-^ + 4A) (3o^ + lOA) X 2 = 20X(Tllm^ 

(10) 

20A(r 1 

(7 + 4 A) (3(7^ + 1 0 A) m 



_20 I ( 7<r Scr 'i 

” 58 m 17^ + 4A 3o^ + lOAJ 



10/Jsina^ sinlSt]^ 

29ml a )’ 

(11) 

where 

j8^ - Y-X 

(12) 


We notice that the mode of speed J(2X) is excited, but does not need 
to be evaluated because it does not affect the middle particle, 

3.6. Radioactive Disintegration of Uranium. The uranium family of 
elements are such that an atom of any one of them, except the last, is 
capable of breaking up into an atom of the next and an atom of helium*. 
The helium atom undergoes no further change. The number of atoms 
of any element in a given specimen that break up in a short interval of 
time is proportional to the time interval and to the number of atoms of 
that element present. If then % X 2 , -^-Xn are the numbers of atoms of 

the various elements present at time t, they will satisfy the differential 
equations 

* We neglect |3-ray products, for reasons that will appear later. 
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du 

dxi 

~dt ~ ~ *^i^u 

dx^ 

— KiXi — K^X^^ 


dXn 
dt 

Suppose that initially only uranium is present. Thus when ^ = 0, w = Wo, 
and all the other dependent variables are zero. Then the subsidiary 
equations are 

{(T-^k)u =cr^lo, 

(a- + Kj) Xi = KU, 

(o* + ^2 = j 

(cr 4- Xfi-i = Kfi^2 ^n-2j 

(rx„ = f(n^lXn-i^ 

The operational solutions are 

(TUo t<0’Uo 

U — ♦ X\ — ^ \ > \ ) 

<r+K ((7- + K)(a + K,) 

KK, CTMo 

“ (o- + /c) (o- + Ki) (o- + /<;) ’ (3) 


(o’ + «)(<^+''l)"- (o’ +««-!)' J 

These are directly adapted for interpretation by the partial-fraction 
rule. In fact 
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Of all the decay constants k is much the smallest. If the time elapsed is 
long enough for all the exponential functions except e-”* to have become 
insignihcant, these results reduce approximately to 

... Mo(1 -«-'<) (5) 

Ki K2 

With the exception of the last, the quantities of the various elements 
decrease, retaining constant ratios to one another. 

On the other hand, if the time elapsed is so short that unity is still 
a first approximation to all the exponential functions, we can proceed 
by expanding the operators in descending powers of a and interpreting 
term by term. Hence we see that at first will increase in proportion to 
t, 0^2 to and to 

In experimental work an intermediate condition often occurs. Some 
of the exponentials may become insignificant in the time occupied by 
an experiment, while others are still nearly unity. We have 


~ K^CT-^Vr^i + •••] ( 6 ) 

and if Krt is small we can neglect the second and later terms in com- 
parison with the first. Thus in this case 

Xr=-Kr^l<r-^Xr^, (7) 

If Xr^i is of the form we can put 




^ (o- + Kj.) (T^ 


(t'+‘ (T + k^ ff*+' 


( 8 ) 


If Krt is small, we can replace the exponential by unity and confirm 
(7). If it is great, 

ft 1 

/ e dt = — + 0 
Jq 


and on continuing the integrations 


KrS\' 


Thus = = (9) 

Kr Kr 

Classifying elements into long-lived and short-lived according as Krt is 
small or large for them, we find that the quantity of the first long-lived 
element after uranium is proportional to ty the next to t\ and so on. 
All i^-ray products are short-lived when t has ordinary values. 

Radium is the third degeneration product of uranium. In rock 
specimens the time elapsed since formation is usually such that the 
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relations (5) have become established. As a matter of observation the 
numbers of atoms of radium and uranium are found to be in the con- 
stant ratio 3*58 x 10“’. This determines Also the rate of break-up 
of radium is known directly : in fact 

1/k^= 2280 years. 

Hence 1 //c = 6*37 x lO^ years. 

This gives the rate of disintegration of uranium itself. 

A number of specimens of uranium compounds were carefully freed 
from radium by Soddy, and then kept for ten years. It was found tliat 
new radium was formed ; the amount found varied as the square of the 
time. This would suggest that of the two elements between uranium 
and radium in the series one was long-lived and the other short-lived. 
Actually, however, it is known independently that both are long-lived. 
The first, however, is chemically inseparable from ordinary uranium, 
and therefore was present in the original specimens ; initially, instead 
of a^i = 0, we have 

K 

For the next element, ionium, we have 

cc^ = Kicr^^Xi = KUQty 

the variation of cci being inappreciable in the time involved. Also 

^3 — K20'~^ 0:2 — ^KK2Uot“. 

Soddy* found that 3 kilograms of uranium in 10*15 years gave 
202 X 10”^^ gm. of radium. Hence, allowing for the difference of atomic 
weights, 

^^3/^0= 7*1 X 10“'^ 

and K 2 = 8*64 x 1 /k 2= 1*16 x 10® years. 

This gives the rate of degeneration of ionium. Soddy gets a slightly 
lower value for l/xg from more numerous data. 

3.7. Some dynamical applications. Suppose we have a dynamical 
system specified by equations such as those in 1.72, save that x is 
replaced by and that the system is at first at rest and then disturbed 
by a force Sr applied to the coordinate The subsidiary equations 
are 

~ 0 (tW 4= r) J ^aCjfiayt^ Sr (?W = (l) 


Phil. Mag. (6) 38, 1919, 483-488. 
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Writing A for the determinant formed by the e’s, and Era for the minor 
of era in this determinant, we have the operational solution 

( 2 ) 

If the determinant A is symmetrical, so that 

Era=E,r ( 3 ) 

we see that a given force Sr applied to the coordinate yr will produce 
precisely the same variation in as the same force would produce in 
yr if it was applied to y^. Thus we have a reciprocity theorem applic- 
able to all non-gyroscopic systems ; linear damping does not invalidate 
the argument if the terms introduced by friction contribute only to 
the leading diagonal of A. 

If the forces reduce to an impulse, so that Sr can be replaced by a-/,., 
the solution becomes 

( 4 ) 

We can evaluate the initial velocities by expanding in descending 
powers of <t. The first term is 

( 5 ) 

where A is the determinant formed by the a’s, and Ara the minor of 
Qrs in it. Hence the initial velocities, found by operating on this with 
<r, are 

= ( 6 ) 


Now the constants Ura are merely twice the coefficients in the kinetic 
energy, which is a quadratic form. Hence the determinant A is S)rm- 
metrical whether the system is gyroscopic or not, and the reciprocity 
theorem for impulses and the velocities produced by them is proved. 

The subsequent motion can be investigated by interpreting according 
to the partial-fraction rule. But let us consider the simple case where 
the system is non-gyroscopic and frictionless. Then 

^ = AU(a^ + a?l ( 7 ) 

where the a’s are the speeds of the normal modes. Then 


ys = 2a 




- Jr sin aty 


where Era{-o,^) denotes the result of putting for a® in En^ The 
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contribution of the a mode to the initial rate of change of is therefore 



An immediate consequence of the presence of the factor Era (- ctO 
the numerator is that if an impulse is applied at a node of any normal 
mode, that mode will be absent from the motion generated. 

Another illustration is provided by Lamb’s discussion* of the waves 
generated in a semi-infinite homogeneous elastic solid by an internal 
disturbance. The normal modes of such a system include a type of 
waves known as Rayleigh waves. These may be of any length, and 
involve both compressional and distortional movement ; if the depth is 
the amplitude of the compressional movement in a given wave is 
proportional to and that of the distortional movement to 
where a and ^ depend only on the wave-length. Lamb found that if 
the original disturbance was an expansive one at a depth /, the ampli- 
tude of the motion at the surface contained a factor but if the 
original disturbance was purely distortional, the corresponding ampli- 
tude contained a factor These factors are the same as would occur 
in the compression and distortion respectively at depth /in a Rayleigh 
wave with given amplitude at the surface. 


CHAPTER IV 

WAVE MOTION IN ONE DIMENSION 


4.1. In a large class of physical problems we meet with the differen- 
tial equation 

•( 1 ) 


dy 

df' 


^ 0 ^-^’ 


where t is the time, x the distance from a fixed point or a fixed plane, 
y the independent variable, and c a known velocity. Let us consider the 
solution of this equation first with regard to the transverse vibrations 
of a stretched string. In this case we know that 

c^^Plm, ( 2 ) 

where P is the tension and m the mass per unit length. Write a* for 
djdti aiid p for d/dx. Suppose that at time zero 

y=m-,%-F{y 


Phil. Tram. A, 203, 1-42, 1904. 


( 3 ) 
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where / and F are known functions of x ; that is, we are given the 
initial displacement and velocity of the string at all points of its length. 
Then we are led by our previous rules to consider the subsidiary 


equation 

{(r^-<^p‘)y=(r‘f{x) + (TF{x) 

(4) 

or 

y 

(5) 

But 

,=i( " ) 

<r — cr \<T — cp <T + cpj 

- i 

(6) 


O’ 1 / O' ^ \ 



CT- — 2cp \(T — cp a + cp) 





2c ^ ^ p 

0) 


In (7) the Ijp has been put last in consequence of our rule that opera- 
tions involving negative powers of p must be carried out before those 
involving positive powers. Hence 


= I {/{^ + cO +/(.r - Ct)}, (8) 


,r^— = 4 f F(x)dx 

o-'-c-’;/ ^ ' 2c^ 'Jo 

•I r- rx + Ct r X ^ ct -i 

■i[f. 

1 fX + Ct 

= i F(x)d^ (9) 

J X - Ct 

and therefore 

3' = H/('^ + cO+/(-J^-cO} + i f F{x)dx. ...(10) 

J X ~ Ct 


This is D'Alembert’s well-known solution. 

This cannot, however, be the complete solution. Equation (1) has 
been assumed to hold at all points of the string ; but if any external 
forces act these must be included on the right of the equation of motion. 
In such problems the ends of the string are usually fixed, and reactions 
at the ends are required to maintain this state ; in the complete equa- 
tions of motion these reactions should appear. They are unknown 
functions of t, and therefore necessitate a change in the mode of 
solution. 
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But D'Alembert's solution can be adapted so as to fulfil all the con- 
ditions. We notice that the solution consists of two waves travelling in 
opposite directions with velocity c. If the initial disturbance is confined 
to a region within the string separated by finite intervals from both 
ends, it will take a finite time before D'Alembert’s solution gives a 
displacement at either end ; hence no force is required to maintain the 
boundary conditions, and the solution mil hold accurately until one of 
the waves reaches one end. 

Again, the initial disturbance is specified only for points within the 
length of the string, that is, for values of .r between 0 and /, say ; and 
by the last paragraph D’Alembert’s solution would be complete if the 
length were infinite. If then we consider an infinite string stretching 
from - 00 to 00 , and disturbed initially so that f (.r) and F{pc) are both 
antisymmetrical about both ^ = 0 and x~l,y will vanish ever after at 
these points, so that no force is required to keep them fixed, and 
D’Alembert's solution for such an infinite string will therefore be 
equivalent to that for the actual string from ^ = 0 to ^ Thus we have 
a formal rule for finding the position of the string at any subsequent 
time : consider an infinite string with the actual values of the initial 
velocity and displacement between x-^ and x = /, but with the dis- 
placement and velocity outside this stretch so specified as to be 
antisymmetrical with regard to both ^ = 0 and = /; then D’Alembert’s 
solution for the infinite string will be correct for the actual string for 
the same values of x. 


4.2. We may approach the problem in another way. Taking the same 
subsidiary equation 





( 1 ) 


except that the effects of the initial velocity will not be considered at 
present, let us solve with regard to ^ as if o* were a constant. The 
justification of this procedure is given later in 4.6. The boundary con- 
ditions are that y = 0 when x = 0 and when x-L Using the method of 
variation of parameters, we assume that the solution is 

y = A cosh — + 5 sinh — , (2) 

c c 

where A and B are functions of x subject to 


A' cosh 


(TX ^ . ,<TX . 

— + B sinh — = 0, 


(3) 
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Substituting in the equation (1) we find 


Hence 

A' sinh ~ + F cosh - =- -f{x) 

c c c 

(4) 

A' 

= ^/(^)sinh^;£' = -^/((.)cosh^. 

(5) 

Now y must vanish for all time when ^-0 ; hence A (0) 

= 0. Thus 


A=f-^/(i)dnh'^di 

(6) 

Also y must vanish when x-L Hence 



iBCOsinh- — J(0cosh- 

(7) 

and 

B(l) = -coth— f -/ (^)sinh-^ d^, ... 
C Jo c ^ c 

(8) 

giving 

B(x) = -coth^ ( — /(^) sinh — 
c Jo c ^ c 


In all we find 



(9) 


y = ( -/($) sinh — (^cosh — - coth — sinh di 
Jo c ^ ' c \ c c c/ 

+ ^ f /(^) sinh 

4 


— ("cosh — - coth ~ sinh 
c \ c c c J 


l 


o-^ sinh or (/ — x)lc 
sinh crljc 

sinh crljc 


di 


Jx C ( 


.( 10 ) 


Before carrying out the integration with regard to we can interpret 
the integrand by the partial-fraction rule. Each integral vanishes when 
o* = 0, and therefore no constant term appears in the solution. But 
sinh (rljc vanishes when 


crl 

— = riTT, 

c 


( 11 ) 


where r is any integer, positive or negative. Hence 


riTT . , 7'LTr^ . V 1 — 3: 

<r . , <r|8inh<r(/-^)/c X “T T 

-sinh 5 — 

C C smh crljc r=-<» 


cri'rr , . . I 

■ cosh (nTr) - 
^ ' c 


^nirctll 


i 


2 2 . nri . nrx ritct 


.( 12 ) 
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By symmetry the corresponding factor in the second integral has the 
same value. Hence 


/ (w sin y di . sm ^ cos (13) 

This is the solution given by the method of normal coordinates. Putting 
t zero, we obtain the Fourier sine series 

y (-^) = I ^ / (0 sill sin (14) 

The effect of an initial velocity can be obtained by a similar 
process. 

But in practice the solution by trigonometrical series is not often the 
most convenient form. It usually converges slowly; but what is more 
serious is that its form suggests little about the nature of the actual 
motion beyond the fact that it is periodic in time 2//c. To find the 
actual form of the string at any instant it is necessary to find some way 
of summing the series, which may be rather difficult. A more convenient 
method is often the following*. We have seen that the interpretation of 
any operator valid for positive values of the argument is equivalent to 
an integral in the complex plane, the path of integration being on the 
positive side of the imaginary axis. Then a factor cosech (rljc in the 
operator can be written 

1 ^p~<TlfC 

and when interpreted by Bromwich’s rule this will give rise to a series 
in powers of integrand. But since the real part of k is 

everywhere positive on the path, the series is absolutely convergent, 
and integration term by term is justified. Hence it is legitimate to ex- 
pand the operator in this way and to interpret term by term. Then for 

sinh aljc A \ j 

and for ^ 


sinh (Txlc sinh (r{l- t)lc 
sinh aljc 


^p-2<rxlc^ 


(1 (1 4 . p- 2 irl/c 4 . p- 4 ,<rllc ^ •••( 17 ) 

On multiplying out either expression we obtain a series of negative 
exponentials of the forms with a > ^, or with i>a. 


* Heaviside, Electromagnetic Theory^ 2, 108-114. 
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Then ^/(^) di = (18) 

But according to our rule is zero when ^ - (a - i)/c is negative 

and unity when this quantity is positive. Hence the integral (18) is 
zero unless a-ct lies between ii and ^ 2 , and then it is equal to/(a - ct). 

Similarly -«)/<? di= - J/ ($) (19) 

and = 0 if ^ + a 

= 1 if ^ + a, 

so that (19) is zero unless ct + a lies between and ^ 2 , and then is 
equal to/(c^ + a). 

It follows that at any instant the effect of the initial displacement 

/(^) at the point $ is zero except at a special set of points where one 

or other of the quantities a ± is equal to Since a does not involve 
we see that this way of expressing the solution reduces it to a set of 
waves moving in each direction with velocity c. The first three factors 
in (16) give 

( 20 ) 

for values of a: greater than t The first term gives |/(^) at time 
(x-i)jc, the second -i/(0 at time (*r + ^)/c, the third -i/(^) at 
time - ^ + 2 (/ - x)}/c and the fourth ^/(^) at time {;r + ^ + 2 (/- ^)}/c. 
The first term represents the direct wave from ^ to .r, the second the 
wave reflected at ^ = 0, the third that reflected at ^ and the fourth 
one reflected first at a: = 0 and then again at The term in 
the last factor of (16) will give four further and similar pulses, each 
later by 2//c than the corresponding pulse just found. These are pulses 
that have travelled twice more along the string, having been reflected 
once more at each end. Similarly we can proceed to the interpretation 
of later terms as pulses that have undergone still more reflexions. 

The part of the solution arising from values of f greater than x may 
be treated similarly. The interpretation in terms of waves is exactly 
the same. 

4 . 3 . As a particular case let us consider a string of length origin- 
ally drawn aside a distance rj at the point x ~ by so that initially it lies 
in two straight pieces, and then released. Then 

= 0<ir<6| 

ijo = r}(l-‘x)l(l-b) b <x<lj 


( 1 ) 
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.( 2 ) 


.(3) 


and the subsidiary equation is 



The solution that vanishes at the ends is 

y = + A sinh ~ sinh 0 <x<h 

y-i] \ — ^ + A sinh ™ sinh - {I- x) b <x<l 
l~b c ^ 

where the constants have been chosen so as to make y continuous at 
x==b. Also a discontinuity in dyjdx at this point would imply an infinite 
acceleration, which cannot persist. Hence we add the condition that 
dyjdx shall be continuous, which gives 

^7(^4- ^ jcosh ^ sinh ^{l-b)^- sinh ^ cosh ^ i)| = 0 

' (4) 

and on simplifying 

A=-^~ y -T - - r ) cosech — (5) 

(T b {I- b) c ^ ^ 

Thus 


__ rjl (l-b c sinh (rxjc sinh <r (/ - b)lc\ 
^^b{l-b) \ I ^ (T sinh(r//c J 

T}1 


0 <x <b 


]l \b ^ cHin\i(rblcsiin\i(r(l-x)lc] ^ j 
smhai/c / 


..( 6 ) 


y ~ b { l - b ) 

Interpreting by the partial-fraction rule, we notice that the contribu- 
tion from or = 0 just cancels the term independent of o-, while the rest 

2* 


gives’ 


“ 1 . VTrb . riTX met _ , 

2 „ "o sm — 7- sin -J- cos — y- v<x<L (7) 


y - ^ -- I -- ^ I 

Alternatively, we can interpret the solution for 0 < ^ < 6 in exponentials, 
for 

sinh <7.r/c sinh <r(/-6 )/c ^ _ g.2„/c>) 

sinh o-Z/c ^ 


(1 - + ( 


and 

Then 


- = ct t>0. 

(T 


-2vl/c + e-i<rl/c + ( 3 ) 

(9) 


^-<r {b-x)Jc ^ = ct—(b- X)j 


.( 10 ) 


Kayleigh, Theory of Sounds 1, 1894, 186. 
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when this is positive, and otherwise is zero. Thus y retains its initial 
value yixll unaltered until time {h - x)lcy when the wave from the part 
where x>b begins to arrive. After that we shall have 

S' ' b(i- b) tr » - K»« -> + »)} ■■ {(5 - j) *- H"* - *)} ■ 

(H) 

This solution remains correct until ct-h-vxy when the wave reflected 
at ^ = 0 begins to arrive. Thenceforward 

(^g-,(b-x)lc _ ^ = {c# - (i - x)} - {ct - Q) + x)\ = 2x. ...(12) 

Hence y = x-x) = -^-^ (13) 









Fig. 6. 
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The part of the string reached by the first reflected wave is therefore 
parallel to the original position of the part where b<x<L 

When - (6 - ir) + 2 (/ - b\ a wave reflected at = / arrives, and 


This holds until ct = 


7)1 


+ + *)-/} (14) 


b{l-b) 

(b + x) + 2(l-b)-, in the next phase 
b{l-b) \ l) b 


y- 


.(15) 


When ct = 2/, the whole of the string is back in its original position ; 
the term in begins to affect the motion, and the whole 

process repeats itself. We see that at any instant the string is in three 
straight pieces. The two end pieces are parallel to the two portions of 
the string in its initial position, and are at rest. For the middle portion 
the gradient dyjdx is the mean of those for the end portions, and the 

1 7)1/C 

transverse velocity is ± 2 middle portion is always either 

extending or withdrawing at each end with velocity c. [Fig. 6.] 


4 . 4 . A uniform heavy string of length 2/ is fixed at the ends. A 
particle of mass m is attached to the middle of the string. Initially the 
string is straight. An impulse J is given to the particle. Find the sub- 
sequent motion of the particle. (Rayleigh, Theory of Sound, 1, 1894, 
204.) 

The differential equation of the motion of the string is 


df ^ dx^ 


= 0 . 


( 1 ) 


Take ^ zero at the middle of the string. By symmetry we need consider 
only the range of values 0<x<L Suppose that the displacement of the 
particle is rj. When t = 0, y and dy/dt are zero except at The 

subsidiary equation therefore needs no additional terms. The solution 
is therefore 


sinho-(/ — ip)/c 
^ ^ sinh (rljc 


( 2 ) 


If now P be the tension in the string, the equation of motion of the 
particle is 


m 


d^Ti 

W 



^3) 


If p be the mass of the string per unit length, 


P = pc\ 


( 4 ) 
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On account of the initial conditions the subsidiary equation tor the 
particle is 



= - 2Prf - coth —-\-Jcr 

c c 

(5) 

Hence 

Jc 

Yj — 


mae + 2P coth o-//c 



Put 

pi mass of half the string 

m mass of the particle 

(7) 

Then 

2P = kme^jl 

(8) 

and 

_ IJjme 

^ ((rlje) + k coth (<t//c) ’ 

(9) 


To interpret this solution by the partial-fraction rule, we recollect that 
the system is a stable one without dissipation, and therefore all the 
zeros of the denominator are purely imaginary. Putting 


crljc-iisi (10) 

we see that the zeros are the roots of 


0) = ^cot <0 (11) 

There is a root between every two consecutive multiples of tt, positive 
or negative; and the roots occur in pairs, each pair being numerically 
equal but opposite in sign. Then 


IJ 


u>etfl 


( 12 ) 


^ me L(i)c (1 -f A:cosec^ w) (l/c) 
the summation extending over all positive and negative values of w 


2/J 


me w (1 -h k cosec' w) 


sin (inctll), 


.(13) 


the summation being now over positive values of w. 

If a root of (11) is wtt + X, where n is great and X < tt, we have 


(n-rr + X) tan \ = k (14) 

and \=klmr (15) 


approximately. The series (13) converges like sisiircu, or 2^. It 

therefore converges fairly rapidly at the beginning, but more slowly 
later. Also the w’s are incommensurable, and therefore the motion does 
not repeat itself after a finite time ; thus the labour of computation 
would be great. 


J 


4 
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To express the solution in terms of waves, it is convenient to change 
the unit of time to //c, the time a wave takes to traverse half the length 
of the string. Also Jim may be replaced by K Then 

^ <r + k coth <r (o- + k) ~{(T — k) e~^ 

O' + /; \ o- -H ((r 4- ^ 

The first term is zero for negative values of the time; for positive 
values it is equal to 

(17) 

After time 2 the second term no longer vanishes. We have 

_ 1 (T (T 

(V+lO' ” “ T(a'~^k) " 

= ( 18 ) 

and the second term is equivalent to 

{1 (^_2) ^^2 (19) 

The third term is zero for ^ < 4 ; for greater values it is easily found to 
be 

~[\-\l+k{t-i) + k‘^{t- 4)=} (20) 

The process may be extended to determine the motion up to any in- 
stant desired. The entry of a new term into the solution corresponds 
to the arrival of a new wave reflected at the ends. 

4.6. A uniform heavy bar is hanging vertically from one end, and a 
mass m is suddenly attached at the lower end. Find how the tension at 
the upper end varies with the time. (Love, Elasticity, § 283.) 

If X be the distance from the upper end, and y the longitudinal dis- 
placement, y satisfies the differential equation 
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where p is the density, Young's modulus, and jPthe external force 
per unit volume, in this case pg. Put Ejp = and let the displacement 
of a particle before the weight is attached be ^o* Then 

- 4 ^=^ 

When iT = 0, ?/o = 0 ; and if the bar be of length /, di/Jdx = 0 when ^ 

Hence ^ i 7) (3) 

After the weight is attached, we still have 

Sf ^ daf~^^ dx^’ 

while when t-O^y - and dyjdt = 0. Hence the subsidiary equation is 



and the solution that vanishes with x is given by 

A sinhor^/c, (6) 

where A is independent of x. 

If w be the cross section of the bar, the equation of motion of the 
mass m is 

= (7) 

the derivatives being evaluated at ^ The subsidiary equation is 

dy 

mahj = mg + mn^y, - (8) 

which gives on substituting for y from (6) 

(a^ sinh — + cosh — ) A - g (9) 

The tension at the upper end is 

gpl + E(tA jc gpl + (or//c) + mca sinh (cr//c) * 

If h be the ratio of the mass of the weight to that of the bar, 

i- ”,; (11) 

pwl’ Em c' ^ ' 

4-2 
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and the tension is 


gm 

W 


1 


1 


.( 12 ) 


Jc cosh(o'//c) + k {(yljc) sinh crljc J * 

We see that gmlmk is the tension due to the weight of the bar alone, 
and gmlv: is the statical tension due to the added load alone. To evaluate 
the actual tension, we expand the operator in (5) in powers of 
Taking Ijc for the new unit of time, we have 
1 _ 

"IK 

/k(T- 


kd sinh o- 4- cosh a 


{kd + 1) - {kd 


kd + 


r r 
1 . 


1 4"— 
^0*+ 1 




.(13) 


The first term vanishes up to time unity, and afterwards is equal to 

2(1-^--(^-1)/A:) (14) 

This increases steadily up to time 3, when the next term enters. Again, 

/’cr — 1 _ ^ ^ kd 2kd 

(ir +17 ” ~ ^ 

= -\+e-ig‘ + 2{tlk)e-W (15) 

and the first two terms, when < > 3, are equivalent to 

2(1- «-«-!)/*) - 2 [1 - - 2 {(# - 3)//l-} w*] 

= 2e-«-3)/* [1 + (2/yi-) (# - 3) - e-2/*] (16) 

This has a maximum when 


l + g-2/* = ?();_ 3) (17) 

Equation (17) has a root less than 5 if 

4/A:>l + (18) 


which is an equality if ^ = 2*7. Thus for ^ 1 or 2 the maximum tension 
will occur before t = b. If yfc = 1 the maximum tension is when t = 3*568, 
and is equal to 3*266 gml^m — 1*633 times the statical tension. If ^ = 2 
the corresponding results are t = 4*368, 2*520 gmlvi, and 1*680 times the 
statical tension. 

The third term enters at time 5, and afterwards is equal to 

2 [1 - ^ - 5)/4: ^ 2 {(^ - b)lkY j. 

If A = 4, the maximum stress is when t = 6*183, and is equal to 2*29 gmlm. 
The statical tension is 1*25 gm/^, so that the ratio is 1*83. 

Love proceeds by a method of continuation, but the present method 
is much more direct, and probably less troublesome in application. 
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4 . 6 . A general proof that the results given by the operational method, 
when applied to the vibrations of continuous systems, are actually 
correct, has not yet been constructed. It would be necessary to show 
that the solution actually satisfies the differential equation and that it 
gives the correct initial values of the displacement and the velocity at 
all points. The proof that it satisfies the initial conditions would be 
difficult in the most general case. We have seen that in one of the 
simplest problems, that of the uniform string with the ends fixed, the 
verification that the solution is valid for the most general initial dis- 
placement is equivalent to Fourier’s theorem. But for more specific 
problems the operational solution is equivalent to a single integral, and 
the direct verification that it satisfies the initial conditions is usually 
fairly easy. To show that it satisfies the differential equation, it would 
be natural to differentiate under the integral sign and substitute in the 
equation. But in practice it is usually found that the integrand, near 
the ends of the imaginary axis, is only small like some low power of 
1/k (the second in the problems of 4.3, 4.4, and 4.5), and consequently 
the integrals found by differentiating twice under the integral sign do 
not converge. But we can proceed as follows. If a part of our solution is 

y ^ f ) ^ ’ 

and this integral is intelligible for all values of x and t within a certain 
range, we have 

9 ^“ fi 0 h 






2 + 

K 


Also 




O C‘ 


Lim 


^x+ch '^Vx y.r-ch 




— 2 + c dK 

K 


The two integrals are identical before we proceed to the limit, and 
therefore their limiting values are the same. Hence 


and the differential equation is satisfied. This argument does not 
assume that the derivatives of y are expressible as convergent integrals, 
but only that they exist. 
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The argument breaks down at points where the second derivatives do 
not exist; as for instance in 4.3 at the points where the slope of the 
curve formed by the string suddenly changes. At these points there is 
a discontinuity in the transverse component of the tension, so that the 
point has momentarily an infinite acceleration. This is why the velocity 
changes discontinuously when a wave arrives. The momentum of a 
given stretch of the string, however, varies continuously; the difficulty 
is the fault of the representation by a differential equation, not of the 
method of solution. 


CHAPTER V 

CONDUCTION OF HEAT IN ONE DIMENSION 


5 . 1 . The rate of transmission of heat across a surface by conduction 
is equal to -kd Vjdn per unit area, where V is the temperature, k a con- 
stant of the material called the thermal conductivity, and dn an element 
of the normal to the surface. Hence we can show easily that in a uniform 
solid the rate of flow of heat into an element of volume dxdydz is 


kVH^dxdydz, 

But the quantity of heat required to produce a rise of temperature dV 
in unit mass is cdV^ where c is the specific heat, and therefore that 
required to produce a rise dV m unit volume is pcdV, where p is the 
density. Hence the equation of heat conduction is 

l{pcV) = kV^V. ( 1 ) 


If we put 

klpc-h^j ( 2 ) 

is called the thermometric conductivity, and the equation becomes 


dV 

dt 


= V, 


(3) 


In addition, there may be some other source of heat. If this would 
suffice to raise the temperature by P degrees per second if it stayed 
where it was generated, a term P must be added to the right of (3). 

It is usually convenient to write o- for d/dt, and 

= (4) 

The operational solutions are then functions of q ; butg must be expressed 
again in terms of <r before interpreting. 
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5.2. Consider first a uniform rod, with its sides thermally insulated, 
and initially at temperature S, At time zero the end = 0 is cooled to 
temperature zero, and afterwards maintained at that temperature. The 
end iP ~ / is kept at temperature S. Find the variation of temperature at 
other points of the rod. 

The problem being one-dimensional, the equation of heat conduction 


is 


dt 


= 0 , 


( 1 ) 


while at time 0, V is equal to S. Hence the subsidiary equation is 


cFV 




The end conditions are that 

F = 0 when ^ = 0, 

F = /S when x=l.\ 

H,„ce W 

The integrand is an even function of and therefore a single-valued 
function of cr. It has poles where 

ql-± m-rr ; that is, or = - (6) 


where n is any integer. But the negative values of q give the same 
values of o- as the positive values, and therefore when we apply the 
partial-fraction rule we need consider only the positive values. The part 


arising from o- = 0 is 

= 

The general term is 

sinh imr {I — o ')! I 

“ (~A^/iV//^);cosh(m7r)Xf-72A2tW7r) ^ 

sin?^-e-»*»’WP, (8) 

rnr I 

and the complete solution is 

F = Ff + 2 -- sin ~ (9) 

n=inv I J 


If TTht^jl is moderate, this series converges rapidly, and no more con- 
venient solution could be desired. It evidently tends in the limit to the 
steady value Sxjl. 
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But if is small the convergence will be slow. In this case we 
may adopt a form of the expansion method applied to waves*. We can 
write (5) in the form 

+ e-^^+ ...)] ( 10 ) 

For if we interpret this as an integral along the path L, the argument 
of q is between ±^7r at all points of the path, and the series converges 
uniformly. Integration term by term is therefore justifiable, and we may 


interpret term by term. Now 

qx = xa^lhy ( 11 ) 

and by 2.4 (16) 1 - Erf (12) 


F=-SrErf-^+ fl-Erf^^Vfl - Erf ...(13) 

When w is great, 1 - Erf iv is small compared with If then 

is moderate, but large, this .series is rapidly convergent, and can 
in most cases be reduced to its first term. This solution is tlierefore 
convenient in those cases where (9) is not. 


5.3. Om-dimendonal flow of heat in a region infinite in both direc- 
tions. First suppose that at time 0 the distribution of temperature is 
given by 

V^H{x) (1) 

We have just seen that the function expressed in operational form by 


satisfies the differential equation 


-Erf 


2ht^' 


U <1^ 


0 , 


( 2 ) 

( 3 ) 


for positive values of x] and by symmetry it will also satisfy it for 
negative values of x, since the function and its derivatives with regard 
to X are continuous when x = 0. Also when t tends to zero this function 
tends to zero for all positive values of x, and to 2 for all negative 
values. It follows from these facts that the function 

K2-.-«>|[l+Erf^J (4) 

satisfies the dififerential equation for all values of x and all positive 
values of t; and when t tends to zero the function tends to zero for 


* Heaviside, Electromagnetic Theory ^ 2, 69-79, 287-8. 
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negative values of x, and to unity for positive values. Hence this function 
is the solution required. 

Suppose now that the initial distribution of temperature is 

V=Ax\ ( 5 ) 

where /(.r) is any known function. Then this is equivalent to 


r=/" (6) 

•/ — 30 


But the solution for positive values of the time that reduces to if ~ 
when ^ = 0 is 


1 


1 -f Erf 




Hence by the principle of the superposibility of solutions the solution 
of the more general problem is 


y=\l ( 7 ) 

= 1 r % (8) 

2ht- 

Put now |=.r + 2^<^A (9) 

Then V= ^ f f{.x + 2/A) e'^' dX (10) 

sjTr J -V) 

This is the general solution obtained by Fourier. 


5.4. If the temperature is kept constant at .r = 0, but the initial 
temperature is /(.r) for positive values of a\ we may proceed as follows. 
If we consider instead a system infinite in both directions, but with the 
initial temperature specified for negative values of x so that 

/(-.r)=-/(.r), (1) 

we see that the temperature at .r = 0 will be zero at all later instants and 
the solution of this problem will fit the actual one. Hence 

r=: ± r/(i) ji _ _i. r /(^) Ji (2) 

( 3 ) 
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If in particular the initial temperature is everywhere /S, / (^) = 8^ and 

. 1 rx/2ht^ 

V=8.-^ e~^^dk 

=sm-^^ (4) 

2ht^ 


Thus the solution of 5.2 is regenerated. In Kelvin’s solution of the 
problem of the cooling of the earth, 5.3 (10) was adapted to a semi- 
infinite region in this way. 

The value of d Vjdx at the end = 0 is obtainable by differentiating 
the solution valid for a semi-infinite region. In Kelvin’s problem 


and 



( 5 ) 


The same result is found by differentiating (4). 

We notice the curious fact that although the original exact solution 
for a finite region in 5.2 (5) is a single-valued function of or, a square root 
of t appears in the approximate solution for a greatly extended region. 
The reason can be seen by differentiating the exact solution. It gives 




which is a single- valued function of cr. But when we use Bromwich’s 
interpretation we find that 



which is again single valued ; but if is specified to be real and positive 
when K is real and positive, it has a positive real part at all points on 
i, and therefore if I is great coth K^l/h is practically unity. The integral 
is therefore equivalent to 

® 

which is our interpretation of o-^S/ky and is equal to 8lhJ{TTt), We could 
have started by specifying the sign of to be negative when k is positive, 
but then and coth KH/h would both have been simply reversed in sign, 
and the same answer would have been obtained. 


6.5. Imperfect cooling at tlie free end of a one-dimensional region^ 
With the initial conditions of 5.2, let us suppose that the end = / is 
maintained at temperature 8 as before, but that the end = 0 is not 
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effectively cooled to temperature zero. Instead we suppose that it radiates 
away heat at a rate proportional to its temperature. At the same time 
heat is conducted to the end at a rate kdVjdx per unit area. These 
effects must balance if the temperature at the surface is to vary con- 
tinuously, so that instead of having F= 0 at the end as before we shall 
have a relation of the form 

dV 

aF=0atir = 0 (1) 


The equation 5.2 (3) is unaltered, and the operational solution is 

F=>S{1 sinh^ 

where A is to be determined to satisfy (1). Hence 

qA cosh ql-a(l^A sinh ql)^0 

F - / 1 - f ) I 

1 ?cosh(//+asinH/J 


( 2 ) 

( 3 ) 

•( 4 ) 


The roots in o- are real and negative, and we can proceed to an inter- 
pretation by the partial-fraction rule as usual. Or, using the expansion 
in “ waves,” we have 


F= (l - . 


q -h CL 


-2^^ . 




If the length is great enough to make the terms involving inap- 
preciable, we can reduce this to its first two terms, thus : 


V^S 



( 6 ) 


If a is great, the solutions reduce to those of 5.2 ; this is to be expected, 
for (1) then implies that F = 0 when x-Oj which is the boundary con- 
dition adopted in 5.2. If a is small, V reduces to S; the reason is that 
this implies that there is no loss of heat from the end, and therefore the 
temperature does not change anywhere. For intermediate values we may 
proceed as follows. If 


^ + a 


( 7 ) 


Bromwich’s rule gives 




1 [ ak 



K=^X‘. 


( 8 ) 

( 9 ) 


Put 
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The path of integration for X is a curve going from to Re^'^^y 

where B is great, passing on the way within a finite distance of the 
origin on the positive side. Denote this path by M. Then 



1 



..(10) 


1 

TTLj 

UrAt 

..(11) 

But 

1 [ ^ 

- 








-l-Erf-^^ (12) 

2lit- 


The second part of (11) can be written as an integral with regard to /x, 
where 

/X - /\ + all (13) 

along a path obtained by displacing ilT through a distance ah] but the 
integrand is regular between these paths and the route M may still be 
used. The second part is therefore 

~ f exp |/x^^ - /X ^2aht + exp (a^h^t + a^) d/jL 


= - exp + ao’) |l - Erf ^ ^ (14) 

and ^ = 1 - Erf ^ - exp (/ + cuz;) {1 - Erf (^ + y)j (15) 

where ^^xj^ht^] y=^aht^ (16) 

Hence F = /S' [Erf ^ + exp (y-4- cur) {l - Erf (^ + y)}] (17) 

This is the same as Riemann’s solution*. 

The temperature at the end is 

Vx=Q = s exp y^ {1 - Erf y} (18) 


whence the temperature gradient at the end follows by (1). For small 
values of the time the temperature at the free end falls continuously ; 
the temperature gradient there is not instantaneously infinite as in 
5.4 (5). For great values of the time we can use the approximation 2.4 
(21), giving 


F.=o = 


ah J (irt) 




Kiemann-Weber, Partielle Differ entialgleichungerij 2, 95-98, 1912. 
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This is equivalent to one found by Heaviside*. Heaviside gives also an 
expansion in a convergent series, suitable for small values of t, but it 
is probably less convenient than the equivalent expression (18) in finite 
terms. We see from (19) that the longer the time taken the better is 
the approximation to (9F/9.r)a;=o given by the simple theory of 5.2. 

5.6. A long rod is fastened at the end ^ = 0, the other end not being in 
contact with a conductor. Initially it is at temperature 0, but at time 0 
the clamped end is raised to temperature 8 and kept there. Each part of 
the rod loses heat by radiation at a rate proportional to its temperature. 

The differential equation is now 


dV 

dt 


car 


(1) 


where a is a constant. Let us put 


and write the equation 


cr + d^ — K^r^ 


d^-F 

dor 


F = 0. 


( 2 ) 

(3) 


At ^ / there is no conduction out of the rod, and therefore d 
vanishes. Also V= S when ^ 0. The solution is 


yr_ ^^ coshr (/- .r) 
” cosh rl 


(4) 


This can be expanded in powers of The first term, which is the 
only one we require, is 


Put 


F = >SV"^ 



K + = X“. 


(5) 

,( 6 ) 

■(7) 


i L t } a 

* W 

But if \ = a + ju, (9) 

the term in 1/(A. - a) becomes 

2 ^ L “P (r “)) ‘”‘1’ (- t) f 

.J«exp(-f){l-Erfi^^-?} (10) 


* Electromagnetic Theory, 2, 15. 
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as in 5.5 (12). The complete solution is 


r=j4exp(-f){l-Erf*-^} 

H.exp(f)(l-E,f€^}] (1.) 

When a = 0 this reduces to 

K=;Sfl-Erf-^') (12) 

so that the disturbance of temperature spreads along the rod, the time 
needed to produce a given rise of temperature at distance x from the 

end being proportional to x\ If t is small enough to make at^ small, the 
error functions will be practically unity except where x is not great 

compared with 2ht^. For such values of x the exponentials are nearly 
unity. Thus at first the heating proceeds almost as in the absence of 
radiation from the sides of the rod. 

But if at^ is great and xj^ht^ small or moderate, the first error func- 
tion in (11) is practically - 1, and the second + 1. Thus in these con- 
ditions 

V = (13) 

This is seen by a return to the original equation to be the solution 
corresponding to a steady state. This will hold so long as at^ -xl2ht^ 

is large and positive, even if xj^ht^ is itself large. The region where the 
approximation (13) is valid therefore spreads along the rod with velocity 

2ah; if at^ has once become|^large, each point on the rod reaches a nearly 
steady temperature at a time rather greater than xl2ah. 

If at^ is large and xl2ht still greater, the solution is nearly 


V= jl - Erf — 1 (14) 

If further the argument of the error function is large, we can write 
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In the regions that have not reached their steady state, the temperature 
resembles that for the problem without radiation, except that the small 

oc 

factor 1 must be introduced. 

5 . 7 . The cooling of the earth. Cooling in the earth since it first became 
solid has probably not had time to become appreciable except within a 
layer whose thickness is small compared with the radius. It is therefore 
legitimate to neglect the effects of curvature and treat the problem as 
one-dimensional. Radiation from the outer surface must have soon 
reduced the temperature there to that maintained by solar radiation, so 
that we may suppose the surface temperature to be constant and adopt 
it for our zero of temperature. The chief difference from the problem of 
5.2 is that we must allow for the heating effect of radioactivity in the 
outer layers. Suppose first that the quantity P defined in 5. 1 is equal 
to a constant A down to a depth //, and zero below that depth. Take 
the initial temperature to be S + mx^ where m is a constant. Then the 
subsidiary equation is 

(i<x<H\ 

and mx) H<x^ , 

and the solutions are 

A 

V - -TV* + S + mx -f Be^^^ + Ce^^ 
fh^ 

and '^mx-^- De~'i^ 

A term in e"^^ is not required in the solution for great depths, because 
it would imply that the temperature there suddenly dropped a finite 
amount at all depths, however great, in consequence of a disturbance 
near the surface. The conditions to determine B, C, and D are that V 
vanishes at x = 0y and that V and dV/dx are continuous at x = H. 


TT PTiop 

B + C+S+AI<t = 0, (3) 

Ber<ia+ Cefiit+AI<T = De-'i", (4) 

Be-iB-Cei^ = De-iJi (5) 

Solving and substituting in (2), we find 


F=/S^(l-« «*) + + ^ {1 0‘®^sinh 0<x<H 

A 

V = S (I - e~^*) + mx + ~ (cosh qH- 1) «“«* x>H. 



( 2 ) 

H<x. I 
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These solutions involve terms of the form - where a is positive. Tc 
interpret, we write 

f -oa f f ^ K^a\ (Ik 

’h\„ (j ^ ») «p (^■< - t) m 

on integrating by parts. But 

iL“''(*’‘“T)x = ‘'" ® 

and a further integration by parts gives 

i I, (^’» - x) § = S. /„ (“ - &) “P - t) 

= htqe-’^--~r’^. ...(9) 


Xa\ c?\ 






^-aV4/tU, ...(10) 


Thus we can obtain a solution in finite terms, which is easily seen to 
be identical with that previously obtained by a more laborious method*. 

The temperature gradient at the surface is found most easily from 
(6) ; we have 




= 



4 . 

- 2hHqe~^^ - 

= 2h {tlirf e- < - // ( 1 - Erf ( 13 ) 


and therefore 
/W\ 8 


+ + //-//Erf- 


p-/rwn 


This is identical with the result given in earlier works of mine except 
that in these the factor 2 was omitted from the last term. 

Jeffreys, Phil, Mag, 32, 575-691, 1916; The Earth, 84. 
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In the actual problem a considerable simplification arises from the 
fact that H is small compared with 2ht^. On this account we can 
expand the solutions in powers of H, and retain only the earlier terms. 
Thus for the surface temperature gradient we have 






S Allr If 

= j + m+ ( 1 - ] 

h{Trt)^ ^ ^ 2h(irty 

and for the temperature at depths greater than II 


V - S{l — e + mx + ^2 2 • i 

AI1\,/f X AIP 


(16) 


5.71. An alternative possibility is that the radioactive generation of 
heat, instead of being confined entirely to a uniform surface layer, may 
decrease exponentially with depth. In this case the subsidiary equation is 

'’aj -fV — ^S^mx) (1) 


at all depths. We already know the part of the solution contributed 
by fi + mx. The remainder is 



..r_A 



1 


But 

Pijr-d^) 

(3) 

and 




f/— a“ “ vy + a q- a/ 



- i - exp (y^ + ax) {1 — Erf (\ + y)}] 

+ h ~ exp (y**^ - our) {1 “ Erf (A - y)}] .... 

.... (4) 

where 

\ = x/2ht^; y = aht^ 

(5) 


Hence W = { exp (/ - ax) - exp (- cur) } 

+ [(1 - Erf X) - ^ exp (y’ + <w) { 1 - Erf (X + y)} 

- J exp (y“ - cur) { 1 - Erf (X - y)}] 

= - Erf X- ^ exp (y“ + our) {1 - Erf (X + y)} 

+ i exp (y= - cur) { 1 - Erf (y - X)}] (6) 
which is the same as the solution obtained by Ingersoll and Zobel*. 

* Mathematical Theory of Heat Conduction^ Ginn, 1913. 


J 


• 5 
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The contribution of radioactivity to the temperature gradient at the 
surface is 

+ ^ 

" ^ ~ 


When y is great, as it actually is, we have 


/0ir' 

V A 

f, 1 \ 

\ 

Lo 



ah J (irt} 


) w 


5.8. The justification of the method is easier in problems of heat 
conduction than in those of the last chapter, because the integrands 

always contain a factor exp • This tends to zero when k tends 


to c±t 00 in such a way that the integrals obtained by differentiating 
under the integral sign always converge, and can therefore be sub- 
stituted in the differential equation directly. But the integrals for the 
temperature are of the form 




and when we substitute in the equation 


dt 




8'^ V 


the integrand vanishes identically. 


CHAPTER VI 

PROBLEMS WITH SPHERICAL OR CYLINDRICAL 
SYMMETRY 

6.1. So far we have treated only problems of wave transmission or 
conduction of heat in one dimension. If our system has spherical 
symmetry, the equation of transmission of sound takes the form 


= 0 ( 1 ) 

where ^ is the velocity potential, and 

T72* 1 3 , 2 9 ^ 

^ “ r" 8r V 0r/ 0r" r 0r 

( 2 ) 
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The differential equation is therefore equivalent to 

= 0 (3) 

This is of the same form as the equation of transmission of sound in 
one dimension, taking the place of ^ ; but differences of treatment 
arise from differences in the boundary conditions. A similar transforma- 
tion can of course be applied to the equation of heat conduction. 

If the symmetry is cylindrical, we take Cartesian coordinates y, Zy 
and put 

if ^ y = x tan (4) 

and if ^ is independent of and we have 

= l (5) 

tar to \ to/ 

which is capable of no simple transformation analogous to that just 
given for the case of spherical symmetry. This fact gives rise to striking 
differences between the phenomena of wave motion in two and three 
dimensions. 

6.2. Consider first a spherical region of high pressure, surrounded 
by an infinitely extended region of uniform pressure; the boundary 
between them is solid, and the whole is at rest. Suddenly the boundary 
is annihilated. Find the subsequent motion. The problem is that of an 
explosion wave. We suppose the motion small enough to permit the 
neglect of squares of the displacements. At all points 

« 

Initially there is no motion, so that ^ is constant, and may be taken 
as zero. The pressure P is - pd^jdtj where p is the density; this is a 
positive constant when r < a, but zero wlien r > a. Then we can take 
the subsidiary equations to be 

•■'“l (2) 

= 0 r>a] 

The pressure must remain finite at the centre. Hence 

V . CTT 

H> = — + Asinh— r<a (3) 

<r c ^ ^ 

= B exp (- nrjc) r> a (4) 

There can be no term in exp (ar/c) in (4) since it would correspond to 

a wave travelling inwards. Now the pressure and the radial velocity 
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must be continuous at r = a; hence ^ and 0^/9r must be continuous. 
These give 

+ A sinh “ = 6xp , (5) 

-- + -A cosh = B exp (- — )i (b) 

d c c c ^\cj 

whence 

A = ^(c-h aa) e~ ] B = “ ~ “ 20 -^ ^ (7 ) 

Thus outside the original sphere 

H> = — (c - acr) e-^ (c + ao-) (»-+a)/^. . . .( 8 ) 


The associated pressure change is 


II 

1 

- - (r-o)/c _ + 05^ g-<r (.-+0)/^^ 

...( 9 ) 

But ~ ~ oj = 0 when t is negative [ 

= ct-a when t is positive, j 

....(10) 

Hence (““"^1 

1 ^-a{r-a)ic _ Q when ct<r- a 


\<r / 

and 

1 

1 

I 

ii 



- ct — r when ct> r — a 

...(11) 

Similarly r~ + a^ 

1 ^-iT{r+a)ic - 0 when ct<r + a 


and 

- ct — r when ct>r + a 

....(12) 


Hence the pressure disturbance is zero up to time (r-a)/c, when the 
first wave from the compressed region arrives, and after (r + a)/c, when 

the last wave passes. At intermediate times it is equal to ^^{r-ct). 

Thus it is equal to pa/27’ when the first wave comes, - pa/2r when the 
last leaves, and varies linearly with the time in between. The compres- 
sion in front of the shock is associated with an equal rarefaction in the 
rear*. 

Within the sphere the pressure is 


per (“r 
r 


cr" 


(TV I 

{c + acr) sinh — \ 

= P j^l - ^ (^ + . ...(13) 


* Cf. Stokes, Phil. Mag. 34, 1849, 62. 
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This is equal to p up to time (a-r)lc, then drops suddenly to 
p(l-al2r)y decreases linearly with the time till it reaches -paj^r at 
time {a + r)ICj and then rises suddenly to zero. The infinity in the 
pressure at the centre is only instantaneous, for the time the disturbance 
lasts at a given place is 2/-/c, which vanishes at the centre. It is due to 
the simultaneous arrival of the waves from all points on the surface of 
the sphere; at other points the waves from different parts of the surface 
arrive at different times, giving a finite disturbance of pressure over a 
finite interval. Hr <\a, the pressure becomes negative immediately on 
the arrival of the disturbance. 

The behaviour of the velocity at distant points is similar to that of 
the pressure. If u is the radial velocity, 

u ~ d^/dr 

~ ~ — ^jr (14) 

If r is great the first term is proportional to 1/r, the second to 1/rl 
The first is therefore the more important. But the first term is simply 
a multiple of the pressure. Hence there is no motion at a point until 
time (r-a)/c, when the matter suddenly begins to move out with 
velocity a/2r. This velocity decreases linearly until time (r + a)/c, when 
it is -a/2r, and then suddenly ceases. The total outward displacement 
contributed is zero. The second term, however, gives a small velocity 
which vanishes at the beginning and end of the shock, and reaches a 
positive maximum at time r/c. It produces a residual displacement, of 
order ajr times the greatest given by the first term ; this represents 
the fact that the matter originally compressed expands till it reaches 
normal pressure, and the surrounding matter moves outwards to make 
room for it. 

6 . 21 . Consider next the analogous problem with cylindrical sym- 
metry. With analogous initial conditions, the subsidiary equation is 

18/ 8^\ cr 

8'sr \ 8t<r/ 

= 0 'aj>a. 

The solutions are Bessel functions of imaginary argument and order 
zero, and Ko^ct'jtIc), The latter is inadmissible within the 

cylinder, because it is infinite when id- = 0. The former cannot occur 
outside it, for the following reason. The interpretation is to be an 
integral along a route through values of the variable with positive real 
parts, and when 'm is great the asymptotic expansion of Iq {kwIc) con- 
tains exp (fcizr/c) as a factor. Hence the solution would involve exp (ortcr/c) 




70 PROBLEMS WITH SPHERICAL OR CYLINDRICAL SYMMETRY 


and therefore a disturbance travelling inwards. The solutions are therefore 


^ = 0<'m<a 

= BKo a<w. 

Also and d^jdm must be continuous at tzr = a. Hence 


AlJ^] 

\---bkJ~), 

(3) 

\ cj 

0 - \c J 


) 

(4) 

Also we have the identity 




-I,{z)K:{z)-\Iz 

(5) 


Hence we find for points outside 'm = a 

® 

But* Io{z) = - I exp (z cos 0) dO, (7) 

’T Jo 

Kq{z)^ I exp (-Z cosh v)dv, (8) 

Jo 

whence 


a f r ( 

^ejJoJo 


^ + -COS0-- cosh'y ) cosOd^dOdv, 
c c / 


Performing first the integration with regard to k, we obtain a function 
of the form II{t-h), Thus 

^ = -~ f f cosOdOdv (10) 

TTC Jo Jo 

where the range of integration is restricted by the condition that 

ct + acosO-m cosh v>0 (11) 

It follows at once that there is no movement at a place until time 
(tjr - a)lc. Integrating next with regard to v, we see that the admissible 
values range from 0 to cosh~^{(c^ + acos ^)/tcr}, provided the quantity 
in the parentheses is greater than unity. Hence 

( 12 ) 

* The notation is that of Watson’s Bessel Functions, 
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So long as ct>w-a^ there will be a finite range of integration with 
regard to the lower limit is then always zero. If + the 
inequality is satisfied for all values of 0 up to tt, and therefore tt is the 
upper limit. If <'nr + a, the inequality is not satisfied when 0 = tt, and 
the upper limit is cos~^ (is: - c^)/a. The disturbance at any point may 
therefore be divided into three stages, the first until the 

second from then till ct~'m a, and the third later. 

We are concerned chiefly with the pressure. This remains constant 
until = -s: - a ; in the next stage it is equal to 

pa fcoB'~H^-ct)ia cos 6 dO , , 

P= . (13) 

{{ct + a cos 

and in the last to a similar integral with the upper limit replaced by tt. 
By applying the transformation 

ct + a"'sy=-2b; ct + a(ios0 -'m= 2b cos^ (14) 

and integrating on the supposition that b is small, we find that soon 
after the arrival of the wave 


P = hp 




(15) 


When the wave arrives the pressure therefore jumps to and 

then falls by | cja of itself per unit time. The corresponding fraction in 
the spherical problem is c/a. At time 'usjc the pressure is still positive ; 
but when ct=^w^a ox more and 6 is greater than ^tt, the integrand in 
(13) is numerically greater than for the supplementary value of and 
thus P is negative. The passage of the wave of rarefaction is therefore 
indefinitely protracted. To find out how it dies down with the time let 
us suppose that ct is greater than w, and that a is small compared with 
either. Then 


P — hp 


a‘ct 


(16) 


approximately. The residual disturbance falls off like t~^. 


6.22. If the motion was one-dimensional, as for instance if the original 
excess of pressure was confined to a length 2a of a tube, the resulting 
disturbance of pressure would consist of two waves, each with an excess 
of pressure equal to Jp, travelling out in opposite directions with 
velocity c. Comparing the results for the three cases, we see that the 
first disturbance at a given point outside the region originally dis- 
turbed, in each case at the same distance from the nearest point of it, 
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occurs in each case at the same moment. The increase of pressure in 

the one-dimensional problem is in the two-dimensional one 
and in the three-dimensional one In the first case the pressure 

remains constant for time 2a/c, and then drops to zero and remains 
there. In the cylindrical problem it begins to fall instantly, and becomes 
negative in an interval less than 2a/c; it then reaches a negative maxi- 
mum, and dies down again asymptotically to zero. In the spherical one it 
decreases linearly with the time and reaches a negative value equal to 
the original positive one at time 2ajc ; then it suddenly becomes zero 
again. 


6.3. Diverging waves produced by a sphere oscillating radially*. 
Suppose that a sphere of radius a begins at time 0 to oscillate radially 
in period ^rrln. We require the motion of the air outside it. 

The velocity potential ^ satisfies the equation 

( 1 ) 

Initially all is at rest; the solution is therefore 

( 2 ) 


When r = a the outward displacement is, say, ^ sin nt when ^>0, and the 
outward velocity cos nt. Hence 


A 



(r“ + 


(3) 


and 





C^ + 


^.exp {cos«^ sin nt - exp 

c^a? r • {j. 

- r “ “ 1' " + T “ " ” r } 


■exp 


ct-^r-^a 


}] 


(4) 


when ct>r-a. 

The solution has a periodic part with a period equal to that of the 
given disturbance, together with a part dying down with the time at a 
rate independent of but involving the size of the sphere. As there 
is no corresponding term in the problem of 6.2 we may regard it as 


Love, Pfoc. Lond, Math, Soc, (2) 2 1904, 88; Bromwich, ih, (2) 15, 1916, 431. 
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a result of the constraint introduced by the presence of the rigid sphere. 
Its effect on the velocity or the pressure is to that of the second term 
in a ratio comparable with (c/wa)l 

6.4. Aspherical thermometer bulb is initially at a uniform temperature 
equal to that of its surroundings. The temperature of the air decreases 
with height, and the thermometer is carried upwards at such a rate that 
the temperature at the outside of the glass varies linearly with the time. 
Find how the mean temperature of the mercury varies*. 

The temperature within the bulb satisfies the equation 

gJ(rF)-^vr=o, (1) 

where (T~h\f (2) 

That at the outer surface of the glass is Gt, where G is a constant. 
But the glass has only a finite conductivity, so that the surface condition 
at the outside of the mercury is 

® 

A' being another constant. The solution of (1) is 

sinh qr (4) 

RR ( ) gives ^ ( ) 

where a is the inner radius of the glass. 

The mean temperature within the bulb is 

Fo = f r^V dr 


■■ {qa cosh qa - sinh qa) 


_ co sh qa ~ sinh qa . . 

” acrq' Ka sinh qa + qa cosh qa - sinh qa ^ ^ 

In applying the partial-fraction rule, we notice that near (t ~ 0 

y ^ + 

ao-q"^ qa{Ka + ^ K(fa^ + ^q^a^) 

-r ^ ^ 1 ^ 

* Bromwich, PWi. Mag. 37, 1919, 407-410; A. R. McLeod, Phil, Mag. 37, 1919, 134. 
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SO that there is a constant lag in the temperature of the mercury in 
comparison with that of the air. 

The other zeros of the denominator give exponential contributions, 
which are evaluated in Bromwich^s paper. 


6.6. A cylinder of internal radius a can rotate freely about its axis. 
It is filled with viscous liquid, and the whole is rotating as if solid with 
angular velocity wq . The cylinder is instantaneously brought to rest at 
time ^ = 0, and immediately released. Find the angular velocity later, 
(Math. Trip. Schedule B, 1926.) 

The motion is two-dimensional, and there is a stream-function 
satisfying the equation 





( 1 ) 


where v is the kinematic viscosity. Since the motion is symmetrical 
about an axis the right side is identically zero. Put 


- = <r, (T = V7-, = w (m 

at 'uidzsX dm/ 

Initially = 

and the subsidiary equation is 

(Y- - 7^) V“\l/ = - 2r^cDo 

The solution is 


( 2 ) 

.(3) 

•(4) 


il/= AIo (rm) -f BKo{r'si) + (7 log zzr + 2) + (5) 

The velocity must be finite on the axis ; hence B and C are zero. Also 
D is independent of 'uj and therefore cannot affect the motion. 

If /be the moment of inertia of the cylinder per unit length and o> 
its angular velocity, the equation of motion of the cylinder is 

00 ) 

= ( 6 ) 


where is the shearing stress in the fluid. Now 


p, = vp 


?«v 


vpa 


?I3? djf 

I 


dm 


(Pl\ 

\m dm/ 


•(7) 


evaluated on the outer boundary at the point (a, 0). Hence 


P» = vp — r [m/o (ra) - 2/o' (m)]. 


( 8 ) 
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Since the cylinder starts from rest, 

/(TO) = - 27n/p Ar [m /© (ra) - 2 / 0 ' (m)] (9) 

Also ao) must be equal to the velocity of the fluid. Hence 

ao) - Arlfi' (ra) + awo (10) 

Eliminating A we have 

<0 [{Ka- - 2) Iq (ra) + ra Iq (ra)] = (Hq [ra 4 (ra) - 2lo (ra)].. .(1 1) 

where ^irvpaK^L (12) 

The operational solution is therefore 

_ raI,(ra)-2U(ra) 

^ ^ (Ka - 2 ) iT(ra)^ ra 4 (ra) ^ ^ 

But 

4 (ra) 1 + \7'‘^ar + + . . . ; 4 ' (ra) = im + + .... (14) 

The contribution from o- =0 is found to be 


7 > ^ 

a+vA 


(15) 


say. This is the ultimate angular velocity. The other terms arise from 
the zeros of the denominator. If we write tk for r, and substitute for K 
in terms of 0 )^ we find that these satisfy 

|(cao - (Oj) ^ 4- 2| Ji (ka) - ka Jo (ka) - 0. 

The solution is tlien of the form Wj + The coefficients are 

determinable by the usual method. 


CHAPTER VII 

DISPERSION 

In the propagation of sound waves in air and of weaves on strings the 
velocity of travel of waves is independent of the period. In many problems 
this is not the case ; waves on water afford an important example. 

7 . 1 . Consider a layer of incompressible fluid of density p and depth 
H. Take the origin in the undisturbed position of the free surface, the 
axis of z upwards, and those of x and y in the horizontal plane. Put 

a 0 a 02 02 , 

dt ^ ^ dx dy dy^ ^ * 


( 1 ) 
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The velocity potential O satisfies the equation 


02 ^ 


+ 




( 2 ) 


At the bottom the vertical velocity vanishes. At the free surface, so 
long as the motion is only slightly disturbed from rest, we have 




(3) 


where ^ is the elevation of the free surface. Then we must have 


cr COS r{z + fl) ^ 

r sill vH 


(4) 


The pressure just under the free surface is ~Tp7% where Tp is the 
surface tension. 

But by Bernoulli’s equation it is also equal to ~ pcr^ - F{i), 

and F{t) does not affect the motion and can therefore be omitted. We 
have therefore the further surface condition 


Tr)C (5) 

Combining this with (4) we have the differential equation for C 

{(T- -(g- T't) r tan rll} C = 0 (6) 


If the fluid starts from rest with C equal to 4, a* known function of cc and 
y, the subsidiary equation will have a term 0-^4 on the right, and the 
operational solution will be 

^ taiTr// 

= cosh t {{g - 2Y) r tan r//}i .4 (3) 

In the corresponding problem of a uniform string in 4.1 the coefficient 
of t was simply 'pc. 

Suppose first that the original disturbance consists of an infinite 
elevation along the axis of y, with no disturbance of the surface 
anywhere else. Then 4 can be replaced by />, and (8) is equivalent 
to the integral 

C = “ J cosh t {{g - Tk^) k tan kH}^ cIk (9) 

The integrand has an essential singularity wherever kH is an odd 
multiple of L cannot therefore cross the real axis within a finite 
distance of the origin, but becomes two branches extending to + x above 
and below the axis. 
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7,2 • The method of steepest descents. The elevation of the surface is 
thus expressed in terms of integrals of the type 

6'=/ Qx^t\f{z))dz, (1) 

J A 

where/(^) is an analytic function, and t is real, positive, and independent 
of Put, following Debye, 

( 2 ) 

thus expressing it in real and imaginary parts. If the integral is taken 

along an arbitrary path, the integrand will be the product of a variable 
positive factor with one whose absolute value is unity, but which varies 
in argument more and more rapidly the greater t is. There will evidently 
be advantages in choosing the path in such a way that the large values 
of B are concentrated in the shortest possible interval on it. Now if 

(3) 


z- X ■¥ ty. 


we shall have 


?/- 

dif ’ da? dif 


dx^ 


0 . 


.(4) 


It follows that R can never be an absolute maximum. But it can have 
stationary points, where 


dx dy 


(5) 


and we know that these points will also be stationary points of I and 
zeros of/'(c:). These points are usually called the ^saddle points,' or 
sometimes ‘cols.’ Through any saddle point it will in general be possible 
to draw two (sometimes more) curves such that R is constant along 
them. In sectors between these curves R will be alternately greater and 
less than at the saddle point itself. The sectors wliere R is greater may 
be called the ‘hills,’ those where it is less the ‘valleys.’ If our path of 
integration is to be chosen so as to avoid large values of 7?, it must 
avoid the hills, and keep as far as possible to the valleys. If then the 
complex plane is marked out by the lines of R constant through all the 
saddle points, and A and B lie within the same valley, our path must 
never go outside this valley; but if A and B lie in different valleys, 
the passage from one valley to another must take place through a saddle 
point. In the latter case the value of the integral will be much greater 
than in the former, and therefore interest attaches chiefly to the case 
where the limits of the integral lie in different valleys. 

The paths actually chosen are specified rather more narrowly; the 
direction of the path at any point is chosen so that | dRjds j is as great 
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as possible. If is the inclination of the tangent to the path to the axis 
of oc, we have 






dll 


^=cos^p-+sinf^^, 


( 6 ) 


and if this is to be a numerical maximum for variations in i/f 


0 = - sin 



+ cos if/ 


dR 

hj 


dn ds ’ 


(7) 


where dn is an element of length normal to the path, drawn so that 
when ds is in the direction of increasing, dn is in the direction of y 
increasing. Hence / is constant filong the path. Such a path is called 
a ‘line of steepest descent.' There will be one in each valley. In general 
the limits of the integral will not themselves lie on lines of steepest 
descent, but can be joined to them by patlis within the valleys. 

In general lines of R constant through different saddle points will 
not intersect; and there will be only one saddle point on each line of 
steepest descent. For the former event would imply that R has the same 
value at two saddle points, the latter that / has, and either of these 
events will be exceptional. It follows that as we proceed along a line of 
steepest descent R will rarely reach a minimum and then proceed to 
increase again. For if R had a minimum dRjds would be zero; but dl/ds 
is zero by construction, and therefore the point would be another saddle 
point. Lines of steepest descent usually terminate only at singularities 
of / (z) or at infinity. 

The path of integration once chosen, the greater t is the more closely 
the higher values of the integrand will be concentrated about the saddle 
points. Thus we can obtain an approximation, which will be better the 
larger t is, by considering only the parts of the path in these regions. 


In these conditions we can take 

/(^)-/(^o)+K^-^o)V"(Co), (8) 

where Zq is a saddle point. Put 

\r{zo)\=^A; (9) 

Then on a line of steepest descent 

/iz)=/(z,)-^Af^, (10) 

and we can put 

arg(^~5;o) = a, (11) 
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on the side after passing through Zq, Then 

2 r exp exp (- ^Atr^) d (r exp ta) 

Jo 




.( 12 ) 


To get from ^ to 5 it may be necessary to pass through two or more 
saddle points, with probably traverses in the valleys between the lines 
of steepest descent. Then each saddle point will make its contribution 
to the integral. 

The error involved in this approximation arises from the terms of the 
third and higher orders omitted from (8). Its accuracy therefore depends 
on exp(-iM;'^) having become small before exp I/"' ( 2 ^ 0 )!) has 

begun to (lifter appreciably from unity. Hence 


itA 


G 


n 

Oi. 


1-6.^ 


ly 




.(13) 


must be large. In most cases the approximation is asymptotic, and does 
not represent the first term of a convergent series. 


7.3. In problems of wave motion we often have to evaluate integrals 
of the form 

C =- f (k) exp (K.r - yO y*, (1) 

where and y are known functions of k. As a rule <^(k) is an even 
function, and y an odd one. When k is purely imaginary y is also purely 
imaginary. We reciuire the motion for large values of t, and possibly 
also of 

The function (k) usually introduces no difficulty. It does not involve 
.r or fy and therefore when these are large enough it can be treated as 
constant throughout the range where the integrand is appreciable. 

It is usually convenient to replace k by ik and y by ty, and to consider 


the equivalent integral of the form 

^ = f ^ exp t (ko; - yt) dK (2) 

ZTT J-QO 

The saddle points are given by 

^ — y'^ = 0 (3) 


so that a given ratio xjt specifies a set of predominant values of But 
y' is an even function of k, and therefore if is a saddle point, ~ ko will 
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be another, and if the adopted path passes through either it will pass 
through the other. We may take kq real and positive. Also 

r(Ko) = ~tyo"^ (4) 

Thus if yo" is positive k - Kq has argument - tt on the line of steepest 
descent, and the contribution from kq is 

^(l7?^l) (5) 

The contribution from - kq (where y ' is negative) is similarly 

^ (p^)^ exp - I (ko ar - yoO (6) 

and the two together give 
/ 2 

ViiTyo'”^ V ^ ~ 0) 

Similarly if yo" is negative the two saddle points give 

( 2 

1 ) ^ W cos (ko^- yo^ + Itt) (8) 

TT I yo r 1/ 


7.31. These formulae, due to Kelvin, are the fundamental ones of the 
theory of dispersion. Consider first the cosine factor, and suppose a: in- 
creased by and t by St Then K^^-yo^ is increased by 

^0 Sa^ - yo + (iT - yo' 0 ^^0 ) 

the term in Sk^ appearing because ko is defined as a function of .v and t 
by (3). But the coefficient of Sx^is zero by (3). If then t is kept constant, 
I will vary with x with period 27r//Co; and if x is kept constant, C will 
vary with t with period 27r/yo. Hence is the wave-length, and 
27r/yo the period, of the waves passing a given place. A phase occurring 
at a given place and time is reproduced after an interval St at a place 
such that Sx = yo St\K^. Hence is the velocity of travel of individual 
waves. It may be denoted by c, and called the wave-velocity. 

But #co has been defined by the equation 

x-y^t - 0 

so that a given wave-length and period always occur when xjt has a 
particular value; they seem to travel out with velocity yo', which is 
called the group-velocity. It may also be denoted by C. In general the 
wave-velocity and the group-velocity are unequal, so that a given wave 
changes in period and length as it progresses. They are evidently con- 
nected by the relation 


( 9 ) 
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7 . 4 . Returning now to 7.1 (9) we can separate the hyperbolic function 
into two exponentials, which will represent wave-systems travelling out 
in opposite directions. One of them is equivalent to 


1 

^ / exp L{KX-yt)dK^ (1) 

where k tanh k H. (2) 

When K is small, 

c = (3) 

+ ( 4 ) 

When K is great, 

c = (Tk)^; C=%{TKf- (5) 


In all ordinary cases Tjg is insignificant in comparison with IP. Hence 
for some intermediate value of k the group-velocity is a minimum ; it 
tends to infinity for very short waves, and to a finite limit for very long 
ones. Three cases therefore arise. If xjt is less than the minimum group- 
velocity, there will be no saddle-point on the real axis, and the disturb- 
ance will be small*. If it lies between this minimum and two 

(positive) values of k will give saddle-points, and each will contribute 

to the motion. If it is greater than {gll)^, the only saddle-point will 
correspond to a short wave. The disturbance at a given point will there- 
fore be in three stages. In the first, leading up to time xl{gH)^y only 
very short capillary waves will occur. Then long gravity waves will 
arrive, the wave-lengths of those reaching the point diminishing as time 
goes on. Superposed on them are further capillary waves, their length 
increasing with the time. At a certain moment the wave-lengths of the 
two sets become equal. This corresponds to the arrival of the waves 
with the minimum group-velocity. From then on the water is smooth. 

7 . 41 . Two typical cases therefore arise according as the wave-length is 
large or small compared with the one that gives the minimum group- 
velocity. Take first gravity waves, such that xjt is small compared with 

* It can be shown that the saddle-points are so placed that the relevant one 
contributes an exponential with a negative index to the solution. This is almost 
obvious from considerations of energy. 


J 


6 
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{gH)^ but large compared with the minimum group-velocity. In these 
conditions we can write simply 

^ = c=^(glK)^; dCldK = ^i(glK^)i ( 1 ) 

The solution is then 

(2s) '‘(■j) '“>(«•»-■)'•< + }*') w 

*'''™ rlilj' ® 


and the amplitude increases towards the rear of the wave train like 


7.42. Take next the capillary waves, short enough for gravity to be 
neglected. Then 

c = {TKf; G = dCldK = ^-(TlK)^ ( 1 ) 

At a given instant the amplitude is therefore proportional to or to 

The front of the disturbance is therefore composed of a series of 
capillary waves whose amplitude tends to infinity, and the time taken 
for them to arrive is infinitesimal 

This impossible result arises from the form assumed for the original 
displacement. In taking io~p we assumed that unit volume of liquid 
was originally released on unit length of the actual line ^ = 0 in the 
surface. The mean height of this mass of liquid was therefore infinite, 
and its potential energy also infinite. The system being frictionless, this 
energy must be present somewhere in the waves existing at any instant, 
and infinite amplitudes are therefore a natural consequence of the initial 
conditions. If instead we suppose that the same volume of fluid was 
originally raised, but that it was distributed uniformly between .r = ± /, 
its elevation was 1/2/ in this range. Expressed in operational form this 
gives 

4 = + 0 - - /)] (2) 


The appropriate solution can be found from 7.4 (1) by introducing a 
1 sin k/ 

factor ^ or into the integrand. If the solution already 


found makes kJ small, this additional factor will be practically unity, 
and the same solution will hold. Waves whose length is large compared 
with the extent of the original disturbance will therefore not be much 
affected by its finiteness. 
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But if kqI is large we must consider separately the contributions from 
the terms in and The factor Ifd will give an extra IIkqI in all 
the solutions. The result will be that at a given instant the amplitude 

of a gravity wave will vary like ko ” ^ and that of a capillary wave like 

Waves whose length is short compared with that of the original 
disturbance are therefore heavily reduced in amplitude. 

On deep water the minimum group- velocity is 18 cm. /sec., corresponding 
to a wave-length of 4 '6 cm. and a wave- velocity of 28 cm. /sec. If the 
original disturbance has a horizontal extent of 1 cm. or so, only waves 
with lengths under 1 cm. or so will be affected, and the amplitudes of 
both gravity and capillary waves will increase steadily with diminishing 
group-velocity. A wave of large amplitude will therefore bring up the 
rear, and will leave smooth water behind it. This is observable in the 
waves caused by raindrops and other very concentrated disturbances. 
But if the extent of the original disturbance exceeds a few centimetres 
the capillary waves produced will be very small, and the largest 
amplitude will be associated with a wave whose length is comparable 
with the width of the disturbed region. The largest wave produced by 
the splash of a brick, for instance, has a length of the order of a 
foot. 

7.5. Two exceptional cases may arise in the treatment of dispersion, 
which are both illustrated in the present problem. The validity of the 
approximation 7.2 (12) depends on exp{jJ/(c)} being proportional to 
exi)(~| on a line of steepest descent. If f'\z) has varied by a 
considerable fraction of A before this exponential has become small the 
approximation will not be good. This may happen if A is itself small, 
or if there are two saddle-points close together. Instances occur when 
there is a maximum or minimum group-velocity, or if the group-velocity 
tends to a finite limit when k becomes very small. In the former case a 
value of a^jt a little greater than the minimum group-velocity will give 
two slightly different finite values of k^. In the latter may be small 
enough for the proximity of - kq to affect the contributions of both. 

7.51. Since y is an odd function of k, we may suppose that when k is 
small 

y = Co#c-C3k^ + 0(k®) (1) 

and then 7.3 (2) is equivalent to 
1 7* 

^ — I i/f (ic) exp i{Kx-c^Kt + Ci>^t)dK (2) 

27r j - CO 


6-2 
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There are saddle-points where 

•; ® 

If, as in 7.4, if/ (k) is constant or tends to a limit different from zero when 
K tends to zero, this integral is nearly 

1 r® 

^ = -1 ij/ (P) cos {k(x- C ot) +C 2 K^t}dK 

Jo 

= ~ i//(0) (c-i 0 ” ^ f cos (v^ - ?nr) dv (4) 

^ Jo 

where m = (cot - a:)l{c<it)^‘ (5) 


The integral involved here is called an Airy integral*. It is finite and 
positive, but not stationary, when ?w = 0; it has a maximum when 
m = r28, and oscillates for greater values of m with steadily decreasing 
amplitude. For negative values of m it tends asymptotically to zero. 



graph given here is adapted from Airy’s table. The integral can also be expressed in 
terms of Bessel functions of order ± J. See Watson, 188-190. 
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Consequently the disturbance we are considering produces an immediate 
rise of the level of the water at all distances greater than though 
this rise is very small at great distances. The maximum rise of level is 
where x is rather less than so that it has not travelled out from the 
origin so far as a wave with the limiting velocity would. The maximum 
is followed by a series of waves of gradually diminishing length and 
amplitude, merging ultimately into gravity waves of the deep-water type. 


7.52 . In the case where there is a minimum group-velocity for a finite 
wave-length, let us, with a somewhat different notation from that used 
so far, denote the minimum group-velocity by yo' and the corresponding 
values of k and y by kq and yo . Put 


Then 


r = To + + 0 + l-y^ + ... . 


( 1 ) 

( 2 ) 


1 

C = j 'A (^) exp t (ko^ - yoO exp i{k^X- yo Kj t - lyo' t) . . . (3) 

with an analogous contribution from the negative values of k ; hence 
with the same type of approximation as before we shall have 



cos {koX - yo^eos {k^x ~ y^K^ t - 


2 - i /*^ 

= -'A('fo)cos(«„;r-y,#)(^y„'"0 */ cos (v^ - mv) dv ...(4) 

.'o 

where 7n = ^ (5) 

(l-yo 0 

The solution is therefore the product of a cosine and an Airy integral, 
the interval between consecutive zeros being much greater for the 
former than for the latter. In the neighbourhood of a point that has 
travelled out with the minimum group-velocity the waves have the 
corresponding period and wave-length, but their amplitude falls off 
rapidly towards the rear. In front of this point the amplitude increases 
for a while and then oscillates. 

In both these exceptional cases we notice that the amplitude associated 
with the critical velocity falls off only like the inverse cube root of the 
time, whereas in the typical case it falls off like the inverse square 
root. Hence the further the disturbance progresses the more will the 
waves with the critical group- velocities predominate in relation to the 
others. In waves on water, however, these phenomena are often modified 
by the greater damping effect of viscosity on short waves. 
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CHAPTER VIII 

BESSEL FUNCTIONS 

8.1. The Bessel functions of imaginary argument are defined by the 
expansions 

4(2)= 2 (1) 

where (n + r) ! is to be interpreted as r (« + r + 1) if w is not an integer. 
Put 

( 2 ) 

Then /„(2.c4) = a;-i" 2 ‘ - 

^ r=oH (« + >•)! 

"0 „-(n + r) 

= .r-i» 2 

r=0 

= exp (3) 

where ^ denotes djdx. Hence 

jr'"" ex\)p~^ = In (2x^) (4) 

Differentiating and substituting for x we obtain the familiar recurrence 
relation 

|{2’‘4(c)} = 2«/,._,(2) (5) 

The expression on the left of (4) is equivalent to the integral 

( 6 ) 

if n is positive ; and if the path L is replaced by a loop passing 

around the origin and extending to -x , the result holds without restric- 
tion on n. Hence 

/n (*) = (1 2) /^exp (''J' + J) , (7) 

which is equivalent to Schlafli’s form *. Putting also 

h'^z = \, ( 8 ) 

..have ® 


Watson, Bessel FunctionSy 176-6. 



BESSEL FUNCTIONS 


87 


provided the real part of z is positive. If 



A.=/4 + (/a2- 1)^, (11) 

the sign of the root being determined by the fact that X and /x must 
approach + oo together, we find 

= ^ f exp(fiz) 1 (12) 

giving, in operational form, 

/«(^)= - (13) 

This can easily be verified for n = 0 by expanding in negative powers 
of p. 

8 . 2 - The integrand in 8.1 (12) has branch points at ±1, but no 
other singularities in the finite part of the plane. The path can therefore 
be modified to two loops from - x , each passing around one of these 
branch points. Consider first (r). The loop around + 1 gives a con- 
tribution denoted, in Heaviside’s notation*, by ^//o(;^). We have 

//„(.) = >- /■ exp(Mz)^ %, (1) 

and if c is great the exponential is appreciable only when jR (/j.) is near 
unity. Then we can put 

( 2 ) 

and write 

TTL •> 

ri+i..l. + ^--^(hO^+...lrf.', (3) 

TTt •'(-2v) 5L 2! j 

wliere, if the loop be supposed indefinitely narrow, we must take (- 2v)^ 
positive imaginary on the upper side of tlie patli, and negative imaginary 
on the lower. Hence 



* Electromagnetic Theory^ 2, 453. 
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Similarly we may consider a loop around - 1 and define 

^o(z) = -[ exp(fis) (5) 

the positive sign being taken for the root when - 1 </x < + 1. Then we 
put 

= ( 6 ) 

and obtain 



By substituting the loop integrals for ffo and in the differential 
equation for /q, namely, 


^0 Idl, 

dz^ ^zdz 


1 = 0 , 


( 8 ) 


we find easily that they both satisfy it and therefore constitute a pair 
of independent solutions. The expansions (4) and (7) are asymptotic, 
but a fuller discussion is needed before a limit can be assigned to the 
error involved in stopping the expansions at a given term. Their physical 
interest is that if the variable is umjcy where 'm has its usual meaning 
as a cylindrical coordinate and o- denotes 0/3^, contains as a 
factor, and its interpretation by Bromwich's rule will have exp K(t- ^mjc) 
as a factor in the integrand, where the real part of k is positive. It 
therefore represents a diverging wave. Similarly represents a pure 
converging wave. They are also intimately connected with the two 
Hankel functions and If arg increases till == ly^ where y is 
real and positive, the loops must be swung round till they pass to + 1 oo 
instead of - oo . Then 


H. (.) = H, (.,) = (A)\xp . (3,- i.) (i . ± ^ ...) 

( 9 ) 

K, (z) = K, (ly) = exp - ‘ (y + - ^ + •••) 

= -^IW{y) ( 10 ) 


To express 1^ in terms of Hq and Kq, we suppose that the loop 
around + 1 passes to - oo on the upper side of - 1. The integral around 
this loop is by definition. Taking the other loop next, we see that 
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(ji‘ - 1)^ is negative imaginary between the branch points, and the 
integral is 


f exp(fiz) 
^TTL J -1 ( 


djX 


1 




.( 11 ) 


(- i) (1 - p.^)- 2i 

Thus ^ = + (12) 

Comparing with (4) and (7) we seem to have the anomaly of a purely 
real function being equal to the sum of a purely real one and a purely 
imaginary one. The explanation is that while and are purely real 
by definition when is real and positive, we have had to define Ho by 



means of a loop passing - 1 on the positive imaginary side. When p. + 1 

is imaginary {ijf - 1)^ ha.s a real part, and therefore //„ is not purely 
real. If we had defined H by means of a loop passing - 1 on the under 
side we should have had to reverse the sign of i in the equation 

Ho(z) = 2l,(z)-.Ko(z) (13) 

The reason this phenomenon does not show in the a.symptotic expansions 
is that when z is great enough Ko is smaller than any term in the 
asymptotic expansion of Ho, and therefore cannot affect this expansion. 
The ordinary Bessel function (z) can be defined by 


Joiz) = Io(^) (14) 

= i{//„(c^)+.ir„(c^)}, (15) 

= J{//W(^) + ^„P'(^)}, (16) 

= cos {z - ;Jir) + higher terms (17) 

In operational form it is given by 

7o(^)= . 


( 18 ) 



90 


BESSEL FUNCTIONS 


8.3. In the case of functions of other orders, the extra factor in 

8.1 (12) is {fM -!- (/X* ~ 1)-}“”, which is unity for /x = 1 and for /x = - 1. 
Hence we can define an Hn and a iT,,, the latter in real form, and the 
former with a real asymptotic expansion. If this is done the first terms 
of the expansions will be independent of ??, a known result. But when 
we proceed to the ordinary Bessel function the complex factors affect 
the argument of the cosine ; we have indeed Stokes’s approximation 

C0S(2:- JWTT- Jir) (l) 

8.4. Heaviside’s procedure is to define 74 (.r) by the operator 

V ’ ’ 

in the present notation The unsatisfactory feature of these definitions 
is that whether the fundamental inter])retation of an operator is an 
expansion in powers of l/p or a complex integral, the former operator 
means 2/o (a:) and the latter ± 2*4 But Heaviside applies a process 
equivalent to changing the variable from fx to v as above, and then 
expanding in ascending powers of tiiis new variable. The resulting 
operators are equivalent to H.2 (4) and (7) above; the method really 
amounts to evaluating the portions of the integrals that arise from the 
neighbourhoods of the branch points, and therefore gives //o and l\^ 
correctly. But the operators introduced at the start do not represent the 
functions Heaviside finishes with. 

8 . 5 . The present expansions can be applied to the solution of the 
problem of 6.21. The operational solution was 



* Kq (ar), as here defined, in accordance with Heaviside’s practice, differs from the 
definition in 6.21, which is that adopted by Watson, following Macdonald; Watson’s 
Jfo is times Heaviside’s. Watson comments on the fact that the extra factor 
obsonres the relation between and the Hankcl functions, but hesitates to remove 
it because the function is already tabulated. But perhaps tables are of less import, 
ance than analytical convenience; when Bessel functions occur in problems they 
more often than not seem to disappear from the final answer. 
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and the pressure is 




which is zero till time (w ~ a)lc, then jumps to {ajvjf^, and proceeds 
to fall by fc/a of itself per unit time. 

8.6. In much of Heaviside’s work use is made of what he calls the 
generalized exponential function, defined by the series 




where w is a proper fraction. The series formally satisfies the differential 
equation 

dx^-y 

and reduces to the ordinary exiionential series when n is an integer. 
ItvS mathematical peculiarity is that the part of the series corresponding 
to negative values of n is always divergent. As usual (>2 + r)! is to be 
interpreted as r(« + r + 1). This series was recently found by Mr A. E. 
Ingham and myself to possess an asymptotic property. Suppose we start 
>vith a given term, say r = - ?«, where m is a positive integer, and con- 
sider the series 

X - wi + r 

.S'- 2 ,, (3) 

r=o(»-W + '-)! 

X 




1 f e'^- 
2Tt j/: k'»- 


(1k 

(k- 1 ) 


if L is so chosen that i k > 1 at all points of it. Then the integrand has 
a pole at k which makes a contribution to the integral. The 
remainder is equal to 

LI. ^ ® 

taken around a loop surrounding the origin and passing to - x at botli 
ends. This is Cijuivalent to 


1 . / 

- sin {m-n)Tr / e 
Jo 
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and the integral is less than 

r dfi = r (m - + 1) 

Jo 

r (n - 7n) sill (n - 111 )^^' 

Hence | S' - ^ | < - ' - — - 

so that the error caused by starting at a definite term is less than 
the preceding term, and decreases indefinitely when x becomes great. 


8 . 7 . As has been indicated, Heaviside’s work is largely concerned 
with the use of as 3 miptotic series, which he justified mainly by appeal 
to actual computation. Less defence of such series is now necessary; 
pure mathematicians have proved their validity in many cases, and 
even used them themselves. Some of Heaviside’s series, however, are in 
their actual application convergent. Thus, if we consider the asymptotic 
expansion of Kq{x\ namely 



\ Hx 


r-\3“ 

2!(8.-r 




3' {2>i -_l_f 


this is formally divergent as it stands ; but in its ajiplication to prob- 
lems of cylindrical waves the argument is crrujc, and the general term 

(allowing for the factor x~^) is 


w» = (- l)’^ 




=(-ir 


r. 3^.. (2/? - 

+ ... Wtt 


, 1 1. 3... 2n-l 

Jtt 71\{271 + 1 ) 2 " 


/ ct 

\2t<7/ 


The general term in the interpretation of {(r-aijc) is therefore 
, \/2 1 . 3 ... 2 ?^- 1 f ct- 

^ TT ^^!(2^^+l)2^‘ \ 2^27 / 


and the series converges like a binomial expansion provided ct<3'oj. 
Though a divergent series is used in the course of the work, it leads to 
a convergent answer. 


8 . 8 . But the operational method raises numerous points concerning 
the relation between Pure Mathematics and Mathematical Physics in 
general. Chapter I of this work suffices to show that the operational 
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method gives the correct results when the system considered has a finite 
number of degrees of freedom; but as soon as we begin to consider con- 
tinuous systems we find that the solution involves operators not 
definable in terms of definite integrations, and the method of Chapter I 
is no longer available for a justification. Bromwich's introduction of 
complex integration serves three important purposes: it enables us to 
prove in large classes of cases that the answers obtained to problems of 
continuous systems are correct ; it provides a formal rule for interpreting 
operators in general; and it gives the answer, in cases where it involves 
new operators, in a form convenient for direct evaluation by contour 
integration. It seems curious that physics, dealing entirely with real 
variables, should find it necessary to use the complex variable to solve 
its problems in the most convenient way. The use of conformal repre- 
sentation in two-dimensional electrostatics and hydrodynamics is of 
course exceptional, and arises from the fact that the general real solution 
of Laplace’s ecpiation in two dimensions is the real part of a function of 
^ + ty ; for other equations there is no such explanation. The real reason 
is that the solutions of the linear differential equations of physics are 
expressible as linear combinations of analytic functions. In the method 
introduced by Bromwich himself, for instance, k is written for d/ct before 
solution, and then for the general value of k the solution is found that 
satisfies the differential equation and the terminal conditions. Such a 
combination of these solutions as will satisfy the conditions when ^ = 0 
is then constructed in the form of an integral with regard to k; this 
integral is the same as would have been obtained by solving operation- 
ally and interpreting according to Bromwich’s rule. The utility of the 
complex variable then rests on the fact that analytic functions of it 
satisfy Cauchy’s theorem. I do not think, however, that the complex 
integral should be regarded as the definition of the operator. In 2.1 (2) 
the is the fundamental notion, and the integral merely expresses 
a convenient rule for evaluating it. The rules of Chapter I in my opinion 
come nearer to the ultimate meaning of the operators, but their exten- 
sion to the operators that arise in the discussion of continuous systems 
awaits further investigation. 



NOTE 

ON THE NOTATION FOR THE ERROR FUNCTION 
OR PROBABILITY INTEGRAL 


The notation given on p. 26 is one that I adopted in 1916 under the 
impression that it was in general use, and I have since used it in several pub- 
lications *, My definition was recently queried by a correspondent, and I have 
not succeeded in tracing its origin. I have, on the other hand, discovered a 
surprising confusion of other notations. The earliest, due to Gauss f, is 



No name is given to the function by Gauss, and there is no sign that he meant 
the notation to be permanent. Of modern writers, Carslaw, Brunt, and Coolidge 
use this notation. Fourier t gives 

s/rr J R 

but has, so fiir as I have traced, no modern followers. Jahnke and Emde, in 
their tables, use 

2 /*'*• 2 
<P(a:)— - I e"* 

yrr J 0 

and call this the Fehlerintegral, The same notation is widely used by other 
German writers. Whittaker and Robinson § also use it and call the function 
the Error Function. The notation 

rno 

Erf^= 1 dx 

was introduced by J. W. L. Glaisher||, who also used IF 

Erfc.r= j e~^ dx— I Jtt -Ysvix. 

0 

The latter function is also called erf a: by R. Pendlebury^. 

♦ Phil, Mag. 1916, 579-585; 35, 1918, 273; 38, 1919, 718. M.N,n,A,S, 77, 
1916, 95-97. Proc. Roy, Soc, A, 100, 1921, 126-6. The Earth, 1924. 
t Werke, 4, 9. First published 1821. 

:J: Thiorie Analytique de la Chaleur, 1822, 458. 

§ Galcultts of Observations, 1924, 179. 

II Phil, Mag, (4), 42, 1871, 294-302. 

IF Loc, cit, 421-436. 

** Loc, cit, 437-440. 
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Whittaker and Watson* use Qlaisher’s notation with the meanings of Erf 

2 

and Erfc interchanged, while Jeans’s t erf a: is times Glaisher’s Erf. 

Vtt 

Numerous other writers use the integrals, but omit to give any special symbol 
or use only the non-committal “ 1.” 

It can hardly be denied, in view of the wide application of the error function 
to thermal conduction, the theory of errors, statistics, and the dynamical theory 
of gases, that it merits a distinctive notation. It is equally clear that none of 
those yet used has obtained general acceptance. Of them, it seems that those 
involving only the single letters 0, ^ are wholly undesirable. In dynamical 

problems these letters are in continual use for couples, while is often wanted 
for a velocity potential and ^ or yj/ for a stream function ; is of course also 
often needed for an angle. Further, 0 has an established meaning already in 
the theory of elliptic functions, namely, the Theta Function of Jacobi, while 
<I> has another meaning in the theory of numbers. To avoid confusion in some 
of the applications of the function it seems necessiiry to use a combination of 
letters, and the only one in frequent use is Erf. 

As to what function should be denoted by this symbol, the most convenient 
is indicated by the practice of the compilers of tables, such as Jahnke and 

2 fa* 

Emde and Dale, who tabulate — I This is an odd function and 

V^r Jo 

becomes unity when .r = x , two properties that make for analytical convenience 
and are connected with the fact that this form, including the factor 27 r~i, 
usually occurs as such in the solutions of the relevant problems. Accordingly 
I think that, of the various notations proposed, 

2 

Erfvr=“p- I 

\7r J 0 

is the most convenient, and is worthy of wider adoption, though I seem myself 
to have adopted it under a misapprehension originally. 

* Modern Analysis , 1915, 335. 
t Dynamical Theory of Gases^ 1921, 34. 



INTERPRETATIONS OF THE PRINCIPAL OPERATORS 
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references to other papers where problems similar to those of the present work 
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Oliver Heaviside, On Operators in Physical Mathematics^ Proc. Roy. Soc. 
A, 52, 1893, 504-529; 54, 1894, 105-143. Fundamental notions and general 
theory: applications to the exponential function, TayloFs theorem, and 
Bessel functions. The arguments used are in many cases suggestive rather 
than demonstrative, and much in the papers would repay reinvestigation. 

Electromagnetic Theory^ 1, 466 pp., 1893 ; 2, 542 pp., 1899 ; 3, 519pp., 1912. 
Published by “The Electrician,’’ reprinted by Benn, 1922. The original 
edition is the better printed. The applications are mainly to electromag- 
netic waves, but heat conduction is also discussed, and much in the work 
can be extended to waves in general. 

Electrical Papers^ 1, 560 pp., 1892; 2, 587 pp., 1892. 

T. J. PA. Bromwich, Normal Coordinates in Dynamical Systems^ Proc. Bond. 
Math. Soc. (2) 15, 1916, 401-448. Gives the first general justification of the 
operational method for systems with a finite number of degrees of freedom, 
and develops a formally equivalent method applicable to continuous systems. 
A rediscassion of the oscillations of dynamical systems is carried out, and 
presents several advantages over the ordinary method of normal coordinates. 
Several problems in wave propagation are solved. 

Examples of Operational Methods in Mathematical Physics^ Phil. Mag. 
(6) 37, 1919, 407-419. P'inds the temperatures recorded by thermometera 
with spherical and cylindrical bulbs when the temperature outside is 
varying, and solves the problem of the induction balance. 

Symbolical Methods in the Theory of Condxiction of HecUy Proc. Camb. Phil. 
Soc. 20, 1921, 411-427. The principal operators that arise in problems of 
heat conduction are interpreted, and the problem of a sphere cooling by 
radiation from the surface is solved. 

A certain Series of Bessel Functions^ Proc. Bond. Math. Soc. (2) 25, 103- 
114, 1926. Discusses the vibration of a circular membrane with given initial 
conditions. 

Be van B. Baker, An Extension of Heavisidds Operational Method of Solnng 
Differential Equations^ Proc. Edin. Math. Soc. 42, 1924, 95-103. The 
partial-fraction rule is extended to cases where the operand is not H (t). 
See also note on the above paper by T. Kaucky, loc. cit. 43, 1925, 115-116, 
which considers operators in greater generality. 
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on some interpretations. 
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A General Expamion Theorem for the Transient Oscillations of a Con- 
nected System^ Phys. Rev. 10, 1917, 217-225. 


Electric Circuit Theory and the Operational Calculus, Bell Sys. Tech. J. 
4, 1925, 685-761 ; 5, 1926, 50-95, 336-384. A general discussion ah initio. 
The operational solution 


zXt)' 


--A{x) 


is interpreted as the solution of the integral equation 


1 

I A(x)e-<xdx. 
kZ {<) j 0 

This equation is then solved by known rules. Numerous electrical applica- 
tions are given. 

Louis Cohen, Electrical Oscillations in Lines, J. Frank. Inst. 195, 1923, 46- 
58. Altci'nating Current Cable Telegraphy, loc. cit. 165-182. Applications 
of Heavisidds Expansion Theorem, loc. cit. 319-326. Contents sufficiently 
indicated by titles. 

Harold Jeffreys, On Compressional Waves in Two Superposed Layers, Proc. 
Camb. Phil. Soc. 23, 1926, 472-481. Discusses diffraction of an explosion 
wave at a plane boundary, with a seismological application. 

J. J. Smith, The solution of Differential Equatiom by a method similar to 
Eeaviddds, J. Frank, Inst. 195, 1923, 815-850. General theory, with appli- 
cations to electricity and heat conduction. 

An analogy between Pure Mathematics and the Operational methods of 
Heaviside by means of the theory of H-Fumtions, J. Frank. Inst. 200, 1925, 
519-536, 635-672, 775-814. Mainly theoretical, bearing, I think, more on 
the relation of the theory of functions of a real variable to mathematical 
physics in general than to Heaviside’s methods in particular. The ideas are 
interesting and useful, though I arn not in complete agreement with them. 
The last paper contains several physical applications. 

Norbert Wiener, The Operational Calculus, Math. Ann. 95, 1925, 557-584. 
A critical discussion, beginning with a generalized Fourier integral. In 
some cases the interpretations of the operators differ from those of other 
workers. 


Most of the above investigators give Bromwich’s interpretation only a 
passing mention, or none at all, but it seems to me at least as general and as 
demonstrative as any other, and more convenient in application. 
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The following appeared while this tract was in the press : 

H. W. March, Tim Heaviside Operational Calculus^ Bull. Am. Math. Soc. 33, 
1927, 311-318. Proves that Bromwich’s integral is the solution of Carson’s 
integral equation, and derives several rules for interpretation from it. The 
author refers to a paper by K. W. Wagner, Archiv fiir Elektrotechnik, 4, 
1916, 159-193, who seems to have obtained some of Bromwich’s results 
independently. I have not seen the latter paper. 

The following are in the press (July 1927): 

Harold Jeffreys, Wave Propagation in Strings with Continuous and Concen- 
trated Loads^ Proc. Camb. Phil. Soc. Some of the results for continuous 
strings are obtained as the limits of those for light strings loaded regularly 
with heavy particles. It is found that the operator arises as 


Lim 





when I tends to zero through such values as make xjl integral. This gives 
the rule 1.8 (5) in terms of definite integration. 

The Earth^s Thermal History, Gerlands Beitrage z. Geophysik. The 
problem of the cooling of the earth is rediscussed, with allowance for 
variation of conductivity with depth. 
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